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abstract 


A  combination  of  field  and  laboratory  experiments  are  made  in  order  to  expand  mu 
knowledge  of  naturally  occurring  oscillators  boundary  layers.  Chapter  1  describes  field 
observations  of  the  development  ol  wave  driven  boundary  layers  at  the  fluid  sediment 
interface,  lender  the  crest  of  the  wave,  this  development  can  be  idealized  as  an 
identifiable  sequence  of  three  parts.  The  latter  parts  of  this  development  are  ne\er 
observed  to  occur  under  the  trough  of  the  wave  despite  similarities  in  wave  orbital 
velocity  and  acceleration.  It  is  proposed  that  wave  induced  boundary  ventilation,  the 
oscillatory  flow  through  the  surface  of  a  permeable  bed,  may  be  responsible  for  this 
apparent  developmental  asymmetry. 

In  chapter  II,  a  laboratory  study  is  presented  of  ventilated  oscillators  boundary 
layers.  These  are  boundary  layers  arising  from  a  flow  which  oscillates  parallel  to  a 
permeable  bed  which  is  subject  to  oscillating  percolation  of  the  same  frequency  as  the  bed 
parallel  flow.  Measurements  of  boundary  layer  velocities,  bed  stress  and  turbulent  flow 
properties  are  presented.  It  is  observed  that  suction  (flow  into  the  bed)  enhances  the  near 
bed  velocities  and  bed  stress  while  injection  (flow  out  of  the  bed)  leads  to  a  reduction  in 
these  quantities.  As  the  ventilated  oscillatory  boundary  layer  experiences  both  these 
phenomena  in  one  full  cycle,  the  result  is  a  net  stress  and  a  net  boundary  layer  velocity 
in  an  otherwise  symmetric  flow.  While  production  of  turbulence  attributable  to  injection 
is  enhanced,  the  finite  time  required  for  this  to  occur  leads  to  greater  vertically  averaged 
turbulence  in  the  suction  half  cycle.  Turbulence  generated  in  the  suction  half  cycle  is 
maintained  in  a  compact  layer  much  closer  to  the  bed.  These  effects  appear  to  hold  for 
Re  ranging  from  10'  to  10'’  and  for  oscillations  other  than  sinusoidal. 

Field  measurements  of  wave  driven  ventilation  are  presented  in  chapter  III.  The 
measurements  are  compared  against  linear  theory  and  ventilation  is  found  to  be  present 
and  predictable  if  the  local  effective  thickness  of  sediment  is  known.  Ventilation  is  then 
investigated  as  a  possible  mechanism  for  the  maintenance  of  beach  slope.  It  is  shown  that 
when  compared  against  net  bed  stresses  resulting  from  nonlinear  wave  forms  or  bottom 
wind,  the  ventilation  induced  stress  asymmetry  is  an  important  mechanism  for  medium 
to  course  grained  sediments. 
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Fie'd  Observations  of  the  Fluid-Granular  Boundary  Layer 
Under  Near-Breaking  Waves 

Daniel  c.  Conley  and  Douglas  l.  Inman 

Cer 'ct  for  Coastal  Stoics.  Scnpps  Institution  of  Oceanography.  La  Jolla. ,  California 

Numerous  synchronized  um«  senes  from  video  earner is,  pressure  season,  current  meier*.  4na  hoi  fiim 
ineroomeicn  oo  n*rur»<  beaches  sho*  that  boundary  lava  development  under  ihe  crest  of  near -hr faking  »i%t* 
can  be  idealized  as  a  process  composed  of  three  Prunes  regimes  hae  referred  to  as  streaking,  roiling  *od 
plurrung.  The  roiling  and  pluming  regimes  fail  10  develop  under  the  trough.  As  a  consequeoce.  there  u  a 
pronounced  asymmetry  in  instantaneous  sand  transport  and  bouod-ry  layer  phenomena  between  the  w  ave  creat 
and  uough.  However,  laboratory  waves  with  Held  scale  periods  and  w  ave  heights  over  Uun  sand  beds  dc  not 
eth>hii  this  aeat  trough  boundary  layer  asymmetry  ,  indicating  that  a  cniscal  elemrn;  of  similitude  is  absent 
jd  laboratory  experiments.  \ke  suggest  that  a»'t  induced  boundary  veruilauon  is  responsible 


INTRODUCTION 

Over  the  past  several  decades,  considerable  effort  has  bee  a 
devoted  to  tbe  study  of  oscillatory  boundary  layer*  ,n  the 
laboratory  under  widely  varying  conditions  However.  little  uc-rk 
has  been  done  to  study  these  boundary  layers  under  field 
conditions  because  of  logistical  difficulties.  In  this  pa*>er  we 
report  od  field  observations  of  the  development  of  the 
fluid-granular  oscillatory  boundary  layer  We  prc>cot  a 
combination  cf  quantitative  and  qualitative  findings  th3t.  although 
not  fully  understood,  are  recurring  elements  in  the  development  of 
the  fluid -granular  boundary  lay  er  in  the  field.  We  also  attempt  lo 
identify  those  aspects  o'  ’hecreticaJ  and  laboratory  work  that  can 
be  applied  to  these  field  observations 

The  Stokes  solution  for  the  lammar  oscillatory  boundary  layer 
over  a  smooth  impermeable  bed  [Lamb,  1945.  p.  620]  is. 


f 


sin<9) 


■  e  "**ud{0— 
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where  um  is  the  amplitude  of  tbe  free  stream  velocity  just  above 
tbe  boundary  layer,  8  *  (or  -  kx )  is  tbe  phase  of  tbe  free  stream,  0 
is  the  radian  frequency,  k  is  the  wave  number.  5e  *  J 2»  /o  is  the 
Stoker  parameter,  c  is  the  direction  Dormal  to  the  bed.  and  v  is 
kiner.'.'.c  viscosity.  Using  ibis  expression,  tbe  laminar  shear 
stress  can  be  denved.  giving 


sin  ^8  - 


C) 


Notice  that  tbe  shear  stress  leads  the  free  stream  velocity  by  n  /  4 . 
This  solution  has  been  experimentally  venfied  numerous  limes 
(e  g..  SUoih .  1968]. 

As  is  true  of  all  turbulent  flows,  the  problem  of  the  turbulent 


designed  to  measure  me?',  quantities  where  the  mean  a  as  taken  as 
tbe  average  over  a  wave  cycle.  Subsequent  experiments  for  the 
most  part  measured  phase-averaged  quantities  in  order  to  idenufy 
tbe  mean  bebaviui  oicr  a  wave  cycle  and  an^werrd  questions 
relating  to  the  range  of  occurrence  of  turbulence  and  the  effect  of 
roughness. 

Transition  to  turbulence  has  been  predominantly  determined 
through  the  use  of  a  w  ave  fricuon  factor.  fw.  Using  a  form  similar 
ic  that  suggested  by  Putnam  and  Johnson  ( 1 9-4 c ) .  Jcnsson  [196?) 
defined  one  form  of  the  wave  fncuon  factor  as 
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(?) 


where  is  the  amplitude  of  the  bed  stress  It  should  be  noted 
that  this  expression  only  relates  the  maximum  bottom  stress  to  tbe 
maximum  velocity,  regardless  of  when  they  occur  and  provides  no 
information  about  tbe  temporal  behavior  of  the  bottom  stress 
While  this  type  of  formulation  may  he  helpful  for  calculations 
regarding  threshold  values  or  maximal  loading,  it  is  of  limned 
value  when  information  about  the  boundary  layer  behavior 
throughout  tbe  wave  cycle  is  desired  Tbe  approach  used  here  is 
to  identify  a  single  critical  oscillatory  Reynolds  number  fie ,  w  here 
the  fnetion  factors  diverge  from  the  lammar  relation  as  implied  in 
relation  (2).  Tbe  idea  is  that,  like  steady  flow,  flows  with 
Reynolds  numbers  below  that  value  would  be  considered  “laminar" 
and  those  w-iU)  higher  values  would  be  considered  “turbulent  ‘ 
Because  tbe  boundary  layer  thickness  in  oscillatory  flow  is  not 
unambiguously  defined  and  is  difficult  if  not  impossible  to 
measure  in  aD>  situation  but  tbe  laboratory,  wc  use  ihc  oscillatory 
Reynolds  number  fie,  based  od  the  orbital  diameter  d6. 

fie  «  —  W 


vscillatory  boundary  layer  suffers  from  the  classical  closure 
problem,  and  no  analytic  solutioo  exists.  On  tbe  other  hand,  many 
models  for  these  boundary  layers  have  been  developed  (see 
Justesen  (1987)  for  a  review),  and  they  have  been  the  subject  of 
numerous  expcnmeotal  laboratory  investigations.  The  earliest 
experiments  ( Bagnold ,  1946;  Inman  and  Bfrven,  1962]  wxre 
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where  da  is  taken  as  the  value  just  above  the  boundary  layer. 

On  the  basis  of  experiments  and  stability  analysis.  Jcnsson 
[1966]  proposed  fie  «  10*  as  the  value  for  transition  between 
laminar  and  turbulent  flow-s  in  the  oscillatory  boundary  layer  over 
a  sro-K*ih  bed  Id  a  detailed  study  uf  wave  fnetion  factors, 
Kamphuis  (1975}  found  this  value  to  be  the  upper  limit  for  the 
laminar  range  over  a  smooth  bed.  In  a  rcanalysis,  Jonsson  (1980) 
suggests  fie  • *  105  as  the  new  lower  limit  for  turbulence  over 
smooth  beds.  Hino  ex  al.  (198?].  using  a  laser  Doppler 
vclocimetcr  in  air  to  measure  velocity  profiles  through  a  purely 
oscillatory  boundary  layer  over  a  smooth  impermeable  bed  (fie  - 
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Conley  and  Inman:  Field  Observations 

3x10^).  partially  confirmed  this.  They  reported  a  turbulent 
boundary  layer  where  most  turbulence  production  uas  m  the 
decelerating  phase.  They  also  reported  that  the  rurhirirnce  is 
almost  entirely  dissipated  by  the  start  of  flow  acceleration  in  the 
subsequent  half-cycle.  This  raises  questions  about  the  idea  of  an 
oscillatory  boundary  layer  which  *.s  considered  to  be  either  purely 
turbulent  or  purely  laminar  throughout  the  entire  cvcle 

Jensen  e:  cl.  [19S9J  measured  the  boundary  layer  development 
in  an  oscillatory  w  ater  tunnel  over  smooth  ana  rough  beds  where#* 
ranged  from  7. 5x10s  to  6x10*.  They  showed  that  the  boundary 
layer  ever  a  smooth  boitom  could  be  both  laminar  and  Turbulent 
w  ithin  one  half -cycle.  Their  results  show  that  the  iirrut  for  purely 
laminar  boundary'  layers  is  approximately  Re  £  10s  and  the  limit 
for  purely  turbulent  boundary  layers  is  approximately  Re  £  10*. 
This  indicates  that  the  idea  of  a  single  'critical*  Re  is 
misconceived  and  explains  some  of  the  previously  observed 
variations  in  the  range  of  the  "critical"  value  of  Re . 

Another  matter  of  interest  bas  been  the  phase  lag  between  the 
occurrence  of  maximum  bed  stress  and  maximum  free  stream 
velocity  um  In  an  oscillating  water  tunnel  study  of  turbulent  flow- 
over  two-dimeDsional  concrete  ripples.  Jonsson  [1965]  found  a  lag 
of  onr  25*  as  compared  with  the  45'  in  laminar  flow.  Re  for 
these  experiments  was  «  6.0x10*  The  ratio  of  wave  displacement 
amplitude.  o0  -  dcl 2,  to  Nikuradse  roughness  parameter. 
was  124.  where  k,  was  determined  by  fituog  measured  boundary 
lay er  velocity  profiles  in  steady  flow  to  tbe  relation  for  the 
turbulent  velocity  profile  ever  a  rough  bottom. 

( u  /  ut )  *  5.o  log(  30;  /  k, )  Here  *  ( T0  /  p  )w  is  the  friction 
velocity  and  rc  is  the  bed  stress. 

Slecih  { 1970]  measured  the  development  of  tbe  oscillatory 
boundary  layer  over  beds  of  loose  sand  in  a  wave  chacnel.  In 
these  experiments,  the  bed  stress  tc  never  exceeded  tbe  threshold 
for  grain  mouoo.  He  found  that  tbe  measured  phase  for  velocity 
over  beds  of  fine  and  medium  saDd  approximately  followed  tbe 
form  of  relation  (1).  However,  for  coane  sand,  where  oc  /  kr  w  as 
about  50.  be  reported  a  phase  lag  wiuch  approached  zero  near  tbe 
bed  (rather  than  tbe  'o’  predicted  by  (1)  and  tbe  25'  measured  by 
Jonsson  [1963]).  Jensen  el  oi.  [1989]  found  that  in  flow  over  a 
smooth  bed  tbe  pbase  lag  $  between  and  um  is  Re  dependent. 
Tbe  phase  lag  is  45*  below  a  Re  of  4x10*,  it  then  falls  off. 
attaining  a  value  of  10’  at  Re  -  10*.  Tbe  phase  continues  to  fall 
off,  though  weakly,  with  Re  above  that  value.  These  results  along 
with  those  of  Sleath  J 1970)  suggest  that  the  phase  lag  is  a  function 
of  Re  and  /  k,  and  that  foe  the  values  of  these  parameters 
expected  as  waves  near  the  sbore.  tbe  phase  will  approach  zero. 

Of  the  aforementioned  laboratory  studies  \Hino  et  al.,  1983, 
Jensen  el  al ..  1989.  Jonsson,  1963:  Kamphuis .  1975.  Sleaih.  1970] 
only  Jonsson  mentions  asymmetry  in  tbe  oscillatory  boundary 
layer  aDd.  in  general,  results  are  only  presented  for  one  half-cycle. 
Al)  tbe  values  measured  by  Jonsson  |J9633,  including  the  free 
stream  velocity,  exhibited  asymmetry  between  the  tw*o  oscillatory 
half-cycles  and  be  attributed  this  asymmetry  to  tbe  mechanism 
driving  tbe  flow*. 

Transition  and  Turbulence  Over  a  Rat  Plate 

Most  authors  agree  that  tbe  first  instability  to  occur  in  tbe 
laminar  boundary  layer  consists  of  two-dimensional 
ToUmien-Scbbcbting  (T-S)  waves.  These  are  small  disturbances 
that  satisfy  tbe  OrT-Somerfeld  Unearned  disturbance  cquation. 
Wheo  tbe  amplitude  of  these  waves  exceeds  a  threshold  value, 
they  begin  to  develop  ihree-dimcnssonality.  manifest  by  periodic 
fluctuations  in  the  cross  stream  direction  At  about  this  same 
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developmental  ume.  secondary  short  wavelength  instabilities  a  e 
excited  b\  these  background  fluctuations,  and  it  is  the  exponential 
growib  of  these  secondary  instabilities  which  leads  to  turbulence 
[eg.  London!  and  Mollo  Chnsler.sen,  1986,  1988)  In 

steady  flow,  the  development  of  T-S  waves  from  inception  to 
three -dimensionality  occurs  o\er  a  distance  approximate!)  equal 
to  five  T-h  wave! digit:,  while  the  development  of  secondary 
instat'd 'ties  aDd  Kreuldown  is  abrupt  and  occurs  ever  one  7  S 
wavelength. 

Conditions  for  T-S  wave  stability  in  steady  flew  have  been 
calculated  by  Jordinson  [1970].  Using  the  Blasius  vcioc;i>  profile, 
be  caku*?led  the  cri’ical  Reynolds  number  Re„  as  570  where 
Re  «  um  S'/v,  8*  is  the  displacement  thickness,  a  ad  u_  is  ihc  free 
stream  velocity.  Here  cnucal  implies  that  unstable  T-S  waves  may 
grow  m  flows  with  Reynolds  numbers  above  this  '  alue  Jordinson 
calculated  the  parameters  for  tbe  first  unstable  T  S  wave  ic? 
develop  as  aS*  *  C.30.  cS *'um  *  0  )2.  c/um  =040,  where  a  is 
tbe  wave  number,  0  is  the  radian  frequency,  and  c  is  the  phase 
speed  Using  tbe  fact  .bat  for  tbe  Blas.js  profile  the  displacement 
thickness  is  approximately  one-third  tbe  boundary  layer  thickness 
6.  this  gives  as  a  streamwjse  wavelength  lt  *  78  Tbe  initial 
cross-flow  wavelength  of  three-dimeDsiooabry  /;  is  approximately 
equal  to  lt.  Relative  to  tbe  scales  of  turbulent  motion,  these  axe 
clearly  long- wave  phenomena 

Considerable  study  bas  been  invested  in  recent  yeas  to 
determine  the  role  cf  turbulent  bursting  m  geophysical  flows  (see 
KUne  et  al  [1967]  for  a  description  of  turbulent  bursting) 
Investigators  report  observations  of  turbulent  bursting  in  boundary 
iaxers  over  plowed  fields  [Mercerei,  J972],  m  tidal  currents  in  the 
Insh  Sea  \Heaiher:ho*\  1974).  and  m  river  flow  where  its 
signature  is  said  to  be  the  boils  present  on  the  surface  of  the  fiow 
[Jockscn,  1976), 

Additionally,  many  investigators  bave  found  that  different 
characteristics  of  the  bursting  cycle  have  roles  in  different  aspects 
of  sediment  transport.  In  unidirectional  fiow,  Weedman  end 
Slmgerland  [1985]  showed  that  sand  streaks  formed  at  lew  friction 
velocity  (ii_<  4  cm's)  in  fine  to  medium  sand  exhibit  the  same 
lateral  spacing  and  statistical  distribution  as  tbe  low  speed  streaks 
found  in  turbulent  bursting.  The  dimensionless  spacing  s4  is 
sut 
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where  s  is  tbe  mean  dimensional  spacing  of  the  streaks 

Relatively  few  field  studies  of  the  flu  id -granular  oscillatory 
boundary  layer  over  permeable  beds  bave  beer  reported  in  the 
literature.  Inman  [1957]  studied  tbe  formation  of  wave-gcneTaied 
ripples  over  sand  bottoms  Near  tbe  surf  zone,  he  observed  npples 
of  low'  steepness  which  were  "ephemeral"  in  nature,  forming  and 
disappearing  in  the  time  scale  of  individual  waves  Sediment 
motion  during  tbe  presence  of  these  ripples  was  described  as  a 
"dense  layer  of  suspended  particles."  fully  distinct  from  the 
discrete  clouds  of  sediment  shed  from  tbe  crests  of  vortex  ripples 
of  more  permanent  form.  Madsen  [1974]  reports  divers' 
observations  that  the  sediment  bed  seemingly  "explodes"  just  prior 
to  the  passage  of  tbe  crest  of  near-breaking  waves  No  such 
phenomenon  was  observed  uoder  the  trough. 

In  a  field  study  using  a  bottom  profiler.  Dmpler  and  Inman 
|1976]  determined  the  range  of  occurrence  of  these  npple  types  in 
terms  of  the  wave  form  of  the  Shields  numbcT  6  where 
pu  3 

e - — —  <*) 

and  p#  is  the  sediment  density  Their  results  indicate  that  npples 
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of  more  permanent  form,  vortex  npptes,  occur  after  the  onset  of 
motion  up  to  0  «  40  Transition  ripples,  the  ephemeral  apples  of 
In/nan  {(957).  occur  for  a  range  of  Shields  number  of 
40  <0  <240  For  values  of  0  greater  than  240.  the  bed  becomes 
Hat  and  oo  ripples  occur  While  studies  by  BagnoiJ  [  I94hj.  Inman 
and  Bonin  [1962}.  Inman  and  Tunsmll  [1972 j.  and  Tunsial!  and 
Inman  (1975)  describe  quite  welt  the  flu  id -sediment  interactions 
over  vorux  ripples  at  conditions  characterized  by  the  lower 
Shields  number  (0540),  no  field  studies  of  the  rtuid-granuiar 
oscillatory  boundary  layer  development  m  conditions  of  higu 
bottom  stress  over  permeable  beds  have  been  reported 

Experiments 

Experimental  Procedure 

Beginning  in  me  early  1980s.  a  senes  of  experiments  designed 
to  investigate  the  fluid- granular  oscillatory  boundary  laser  m  the 
region  of  high  bottom  stress  (9  *  240  j  were  conducted  adjacent 
to  the  ocean  pier  of  Scripps  Institution  of  Oceanography  m  water 
depths  of  t  -  3  m.  All  the  measurements  were  ai3de  just  seaward 
of  the  breakpoint  of  the  largest  waves  m  a  zone  where  wave 
heights  were  0.5  to  1  m.  and  the  slope  of  the  fine  sand  beach  was 
approximately  1  in  50  The  sand  has  a  median  diameter  O  of 
about  150  um  and  is  predominantly  quartz  wub  atvui  I0rr  by 
weight  of  bwavy  minerals,  mainly  hornblende 

Sensors  were  mounted  on  a  tn angular  “surf  spider"  constructed 
of  welded  aluminum  tubing  with  a  spacing  of  I  85  m  between  legs 
(Figure  1)  The  surf  spider  was  floated  into  place  by  divers  using 


aj  bags  and  was  securely  augured  into  the  bed  by  means  of 
helical  plates  mounted  on  rods  extending  through  the  hollow  legs 
The  spider  was  oriented  with  two  legs  along  a  cross -shore  hnc 
and  the  sensors  were  mounted  u>  avoid  inter Jerence  with  each 
other  and  uni  the  legs  ( Inman  el  al . 

The  sensors  were  mounted  m  a  plane  perpendicular  to  me  cross 
shore  direction  and  included  a  pressure  scosot.  a  two- component 
electromagnetic  iem)  currem  nKWri.  and  two  cylindrical  ot  fur. 
probes  which  were  mounted  with  a  1cm  separation  so  the 
longshore  an1  vertical  directions  Prior  to  each  run.  the  probes 
were  adjusted  so  that  the  bottom  probe  was  approximately  1  err, 
above  the  stiil  bed  level  {Figure  !)  la  addition,  an  8-mm  video 
camera  was  aligned  and  focused  so  that  the  hot  film  probes  and 
the  surrounding  sand  bottom  were  m  the  field  of  view  Numerous 
suii  photographs  were  taken  during  the  experiments 

The  pressure  sensor  wi»  a  temperature  compensated  Statham 
model  PA506-3?  absolute  pressure  transducer  The  cm  currcr! 
meter  was  a  Marsh-McBuney  model  512  OEM  with  a  two axiv 
4-cm-diarreter  spherical  probe  with  a  nearly  constant  response  up 
to  l  0  Hz.  The  current  meter  was  placed  about  50  cm.  above  the 
bed  and  used  to  measure  the  instantaneous  free  stream  wave 
velociry  The  current  meter  was  also  used  :n  in  situ  calibrations 
or  i he  hot  fdms  The  aneraomeuy  was  composed  of  T 5 1 
Incorporated  model  1750  constant  temperature  anemometers 
< electronics}  along  with  TSI  Incorporated  model  »2h'’-b0NuCl 
cyiindncai  hot  film  probes  The  probes  extended  from  a  5  2-mrv. 
diameter  support  rod  and  could  be  raised  to  the  same  level  a>  The 
current  meter  for  signal  comparison  or  lowered  to  within  1  cm  of 
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Fig  I.  The  rurf  jptder  platform  used  for  jervsor  deployment  The  hot  film  package  is  mounted  on  a  rack  and  pinch  roller  f^r 
verucal  adjustmenu  in  the  field 
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f'S  2.  Typical  limf  serif*  from  luce  5  3^89.  experiment  Sensor  outputs  Hop  to  hotionii  *s<  pie.«<u;e  rr^d  d/pih  crr«*  shcr? 
'eloaty  (tm  current  meie?  »nd  bouom  veltxnues  (hot  Dims)  Hot  him  pioPer  1  and  2  were  approximate!!,  0  5  cm  and  3  5  cm 
above  th-  'at  ren*  bed 


the  bed  The  convention  used  designates  onshore  How  as  positive 

The  video  camera  was  a  Sony  V9  S mm  Video  camera  and  w as 
used  id  data  rids  of  approximately  10  min  Tbc  video  was 
synchronized  v*  ith  the  data  through  tbe  use  of  the  time  signal 
recorded  on  tbc  video  record  Prc-  and  posicahbrations  indicated 
that  tbc  video  and  data  were  synchronized  to  within  e  ODe  frame 
(1/30  sj. 

Data  from  all  the  sensors  were  cabled  from  the  surf  spider  up 
onto  the  pier  and  recorded  af  12$  Hz.  A  typical  section  of  data  js 
shown  in  Figure  2.  The  video  and  sull  photos  were  used  to 
determine  boundary  layer  development  patterns,  while  the 
synchronized  data  records  were  used  to  determine  onset  conditions 
and  to  verify  interpretations  of  the  various  patterns  m  the 
development  of  the  fluid-granular  boundary  layer. 

We  shall  rely  primarily  on  the  June  5,  1989,  experiment,  which 
provided  the  most  comprehensive  data  scl,  Tor  most  of  the 
quantitative  results  presented  here.  During  this  experiment 
incident  waves  bad  a  spectral  peak  T  of  12.8  s.  a  Significant  wave 
height  of  70  cm,  and  mean  water  depth  h  of  190  cm.  For 
significant  waves,  this  corresponds  to  a  fie  of  1.2  x  10*.  a  Shields 
number  of  260  and  a  ratio  of  (at  /  kt )  *  5400.  Using  the  results 
of  Jensen  el  al  (1989,  Figure  34  j,  this  corresponds  to  a  boundary 
layer  thickness  of  6  *  3.0- 3.1  cm.  Six  5J2-s  data  runs  were 
collected  over  a  2  hour  period.  Tbe  power  spectra  of  cross-shore 
velocities  for  all  runs  are  shown  m  Figure  3.  Supplementary 
calculations  nave  been  made  from  underwater  photographs  taken 
dunog  two  experiments  on  August  8.  1985,  and  September  6, 
1985,  when  T  *  6.4,  9.4  s;  a  69,  71  cm;  and  h  *  120, 
150  cm,  respectively.  These  correspond  to  Re  =  0.9x10*,  1  1*10*; 
9  *  420.  360 ;o/kt  *  3400,  4500;  and  8  =  1,9-2  0.  2  6-2  7  cm. 


Dc:a  Analysis 

Each  of  tbe  six  data  runs  for  /une  5.  1989  consisted  rf  512  < 
of  data  recorded  at  128  Hz.  Laborator>  calibrations  were  applied 
to  the  etn  current  melt'  and  pressure  sensor  records  Tbe  hot 
films  were  calibrated  following  Flick  end  Geer  ye  {1990]  m  whose 
work  tbe  cm  current  meter  is  used  as  a  velocity  standard  io  field 
calibrations.  From  calibration  runs  in  which  tbe  bot  film  probes 
and  current  meter  were  at  the  same  elevation  in  a  plane  parallel  to 
tbe  wave  froDt,  a  nonlinear  hot  film  calibration  following  King's 
law  was  calculated  by  comparing  the  hot  film  and  current  meter 
signals  Once  cabbrauons  were  applied,  the  hot  film  records  were 
derectified,  rendering  them  of  value  for  computing  mean  flow  s  as 
well  as  instantaneous  speeds.  Spectra  for  each  ran  were  calculated 
by  eosembbng  coefficients  c.  ''ulaied  from  four  data  blocks  of 
128  s. 

The  videos  were  first  observed  in  the  ir  eruirery  in  order  to 
develop  impressions  of  boundary  layer  development  patterns 
When  certain  events  or  sequences  were  identified,  tbe  entire  video 
record  was  systematically  analyzed  using  stop  aclioo  or  frame 
stepping  to  identify  tbe  points  of  occurrence  or  to  measure  specific 
length  scales.  When  times  of  occurrence  were  identified,  cross 
reference  could  be  made  with  tbe  velocity  and  pressure  signals. 
When  events  or  phenomena  were  observed  in  the  videos  but  could 
not  be  resolved  with  "freeze"  frame,  still  photographs  were 
employed.  While  sull  photography  provides  good  quanutaUve 
information  about  length  scales,  event  specific  flow  parameters 
could  not  be  obtained,  as  there  was  do  synchronization  between 
tbe  sull  camera  and  sensors. 

Tbe  natural  variability  of  wave  amplitude  and  period  resulted  in 
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Fig.  ?.  Six  po*er  spectri  with  8  DOF  each  of  oosj  shore  v«)oc>ues 
measured  b>  em  current  meter  on  Jooe  5.  1989. 


wsdely  varying  condiuons  which  significantly  impacted  the  video 
records.  Tbe  most  dramatic  effects  occurred  following  the  passage 
of  groups  of  higber  waves  in  which  ciouds  of  suspended  sediment 
obscured  ibe  bottom  for  several  wave  half -cycles  Other  cycles 
were  rendered  unusable  by  tbe  presence  of  kelp  or  because  the 
wave  intensity  barely  exceeded  tbe  threshold  of  notion.  Tbe 
3072  s  of  record  from  tbe  June  5  experiment  contained  197 
individual  waves  of  sufficient  magnitude  to  set  the  lop  layer  of 
sediment  in  mouon  at  some  point  under  the  crest  of  the  wave 
Tbe  number  of  waves  in  which  the  bottom  was  dearly  visible 
throughout  all  aspects  of  tbe  fluid-granular  boundary  layer 
development  was  less. 

Bound  as  y  La  ye*  Development  Unde*  die  wave  Crest 
Streaking 

Beginning  with  the  fluid  at  rest,  tbe  fust  visual  indication  of  the 
fluid-gjamila/  boundary  layer  development  is  tbe  initiation  of 
motion  where  grains  begm  to  roll  along  the  bed  and  then  rapidly 
develop  into  a  series  of  sand  streaks  (Figure  4).  Although 
occurring  in  apparently  laminar  oscillatory  flow,  the  streak 
development  resembled  the  sequence  described  by  Weedmcn  and 
ShngerUmd  [  1985}  for  tbe  formation  of  streaks  tn  fine  and  medium 
sand  under  conditions  of  weak  friction  velocity  Tbe  sneaks  here 
are  visually  manifested  by  the  dark  lines  which  result  when  quaru 
grains  (large  white)  are  separated  from  hornblende  (small  dark) 
along  lateral  shear  lines  by  a  sorting  mechanism  described  by 
Inman  el  of.  (1966).  These  streaks  initially  appear  as  rather  long 
O(]0  cm)  and  regularly  spaced  lines  which  become  more 
entwined,  giving  an  echelon  or  lambda-like  pattern.  For  tbe 
durauoo  of  this  regime,  tbe  sediment  is  to  a  rolling  gram  or 
shooting  flow  where  some  portior  of  the  immobilized  bed  »* 
visible  and  the  thickness  of  the  layer  of  sediment  in  motion  is 


never  more  iban  a  few  gram  diameters  This  entire  regime  ss 
labeled  'streaking  " 

Visual  impressions  of  the  particle  paths  of  neutrally  buoyant 
organic  material  and  sediment  grains  remaining  sn  the  fluid 
colunm  from  previous  events  suggej*  that  the  fl->w,  starting  from 
rest  after  the  reversal  of  mouon.  is  laminar  As  the  flew 
accelerates,  it  attains  the  magnitude  necessary  for  the  onset  of 
gram  motion,  u  .  and  grams  of  sediment  on  the  surface  cf  the  bed 
begin  to  roll  in  ’he  stream “ise  direction  Tbs  js  quickly  followed 
by  the  appearance  of  sand  streaks  which  can  persist  for  the  entire 
half-cycle  of  wave  mouon  when  the  motion  is  weak 

Analysis  of  two  photographs  taken  at  times  of  “no  morion* 
during  tbe  August  8  experiment  gives  a  mean  spacing  fen  151 
remnant  screaks  as  4.8  mm  with  a  standard  delation  of  1.7  mm 
Tbe  mean  spacing  for  69  streaks  from  a  photograph  taken  during 
conditions  of  weak  motion  on  the  September  6  experiment  w  as 
3.9  mm  with  a  standard  deviation  of  1.3  mm.  Using  the  wave 
parameters  for  August  8  and  the  assumption  that  streaking  occurs 
around  or  ■  15*.  the  results  of  Jensen  el  ai  (1989}  give  a  'fiction 
velocity  of  2  cm/s.  This  results  sn  a  normalized  streak  spacing  r‘ 
of  89  For  September  6,  u.  =  1.2  cm  and  r*  =  43  As  these 
streaks  occur  during  apparently  laminar  flow,  u  is  unlikely  that 
they  are  tbe  turbulent  low  speed  streaks  observed  by  Weedman  and 
Slingerlond  [  1985).  This  is  supported  by  the  fact  that  the  observed 
streak  spacing  in  oar  experiments  is  lower  than  thai  expected  for 
turbulent  streak  spacing  in  unidirectional  flow- 

Roiling 

Tbe  next  regime  begins  with  the  abrupt  Iran:  tion  of  the  fluid- 
granular  boundary  layer  to  a  thicker  carpet  or  sheet  type  flow 
iFigure  4).  During  this  regime,  ibe  layer  of  sediment  in  motion 
resembles  a  carpet  whose  surface  is  composed  cf  tufts  of  sediment 
laden  fluid  The  'tufts*  scxrm  to  be  regular  in  sue  with  heights 
above  the  “at  rest"  bed  of  about  1  6  cm  and  spacing  between  the 
tufts  of  about  6  cm.  Tbe  result  is  a  layer  of  sediment  in  motion 
whose  easily  recognized  surface  is  fully  three-dimensional  w»*h 
variations  id  elevation  in  both  the  cross  and  stream  wise  directions 
(Figure  5j.  Tbe  enure  video  record  from  June  5  was  analy  zed  m 
order  to  isolate  examples  of  this  regime  Of  the  197  “significant" 
waves  io  the  record.  109  half-cycles  undcT  lbe  crest  exhibited  this 
carpet  flow  Of  these,  31  had  an  unobstructed  view  of  the  bottom 
at  the  time  of  transition.  Tbe  abruptness  and  clanty  of  this 
iransiUt.3  is  supported  by  tbe  fact  that  repeat  estimates  of  the  onset 
o,  transition  for  31  different  events  resulted  in  a  mean  variation  of 
onset  ume  of  3  video  frames  (0  1  s). 

Under  relatively  low  amplitude,  long  period  swell  as  sometimes 
prevailed  during  tbe  earber  deployments,  this  carpet  flow  can 
begin  *o  exhibit  stxeamwise  periodicity  much  like  the  wales  in 
corduroy  (Figure  6).  If  the  flow  sequence  terminates  with  this 
flow  regime  and  tbe  flow  has  developed  shreamwise  penodiciry 
the  bed  exhibits  transition  ripples  at  the  cessation  of  motion 
‘(Figure  7).  If.  as  was  tbe  case  throughout  the  iuoe  5  etpciiment. 
the  flow  sequence  terminates  with  this  regime  in  tbe  three- 
dimensional  state  (Figure  5).  tbe  bed  deposited  is  flat.  V*’e  have 
labeled  this  regime  “roiling.* 

Jensen  el  al  (1989]  defined  a  “temporal*  friction  factor  /Jlcsj) 
The  expression  for  this  friction  factor  contained  a  phase  offset  of  v4 
as  required  by  (2).  However,  as  they  showed  in  their  paper,  this 
phase  lag  is  a  function  of  Re  \Jensen  et  al.,  1989.  Figure  1 1 )  so 
we  propose  a  slightly  modified  temporal  friction  factor 
/*(or)  defined  by 
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Fig.  8.  Temporal  friction  factor /  *  versus  Reynolds  cumber  Re  for  coanint  wave  phase  o / .  »  here  90'  correspond*  ic  maximum 
orbital  velocity.  The  heavy  solid  bne  represents  laminar  Stokes  bouoda r>  layer,  and  ihe  heavy  dashed  hoe  u  ihe  Lm»t  for  fully 
turbuleo’  flow.  GrcJes  are  dau  points.  All  are  from  Jcmo t  rt  ai.  { 19S9.  Figure  8]  bones  *rr  po,nu  reproduced  in  o>-?  Figure  9 


Id  tbc  original  expression,  the  phase  lag  0  was  fixed  at  ny 4  but 
here  is  dependent  on  fie  and  equal  to  45'  for  fie  <  4x)04  and  falls 
off  to  *8*  for  Re  =  10'.  Jensen  et  ai  constructed  a  plot 
of  f‘  versus  fie  in  which  turbulent  transiuon  can  be  secD  to  be  a 
fuocuoD  of  both  phase  and  Reynolds  Dumber.  Using  /*  from 
relation  (7),  we  recreate  this  plot  id  Figure  8  Tbe  light  solid  lines 
indicate  bow  the  friction  factor  at  different  wave  phase  (or )  vanes 
from  the  laminar  case  Fully  turbulent  flow  occurs  to  tbe  ngbt  of 
tbe  dashed  line  aod  transitional  flow  is  assumed  to  occur  in  tbe 
region  between  tbe  dashed  and  laminar  line.  Tbe  identified  trend 
is  for  transition  to  occur  at  lower  phase  with  larger  Reynolds 
number. 

We  use  these  results  to  suggest  that  the  transition  from  streaking 
to  roiling  is  the  fluid-gjanular  oscillatory  boundary  layer 
expression  of  tbe  transition  to  turbulence.  This  can  be  seen  in 
Figure  9  which  is  a  plot  of  the  wave  phase  at  which  the  onset  of 
roiling  occurs  versus  tbe  fie  of  the  w  ave  in  which  it  occurred.  In 
order  to  compute  Reynolds  numbers  for  this  plot.  um  was  defined 
as  tbe  peak  velocity  during  tbe  wave  balf-cycle  in  which  roiling 
was  observed.  Tbe  balf-cycle  was  defined  by  subsequent  zero 
crossings  in  tbe  velocity  record,  and  db  was  calculated  by 
integrating  velocity  over  tbe  balf-cycle.  Tbe  wave  phase  was 
calculated  as  or  *  axcsin(  ur  /  um )  where  ut  was  the  free  stream 
velocity  at  onset  of  roiling.  As  can  be  seen,  tbe  transition  to 
roiling  follows  tbe  trend  of  tbc  transition  to  turbulence.  That  is. 
transition  occurs  at  Iowct  phase  for  higher  Reynolds  number.  It 
should  also  be  clear  from  this  figure  that  sneaking,  w  hich  prtctdtt 
roiling,  occurs  at  a  phase  angle  of  15*  or  less.  It  should  be  noted 
that  Figure  8  corresponds  to  flows  over  a  smooth  and  immobile 
bed.  It  is  not  surprising  that  transition  in  a  flow  over  a  loose 
sediment  bed  occurs  at  lower  fie  than  that  expected  under  tbe 
above  idealized  conditions. 

Of  the  88  remaining  waves,  those  with  tbe  most  energetic 
half-cycles  unHi-r  the  crest  were  isolated  for  analysis.  The  highest 
fie  calculate  such  a  flow  was  2.6x10s.  and  29  of  33  haJf- 


cvcles  studied  bad  Reynolds  numbers  below-  2x)(r.  The  lowest 
fie  for  a  wave  jo  which  roiling  was  observed  was  2  2*1 0s  This 
suggests  that  waves  in  which  roiling  is  not  observed  under  the 
crest  cither  experience  transition  late  in  the  balf-cycle  when  ihe 
wave  is  approaching  deceleration  or  not  at  all 


PkaSE  (DEGREES) 


Fig.  9.  Wave  phase  at  the  onset  of  roiling  versus  Ft .  where  a  phase  of 
90’  correspond*  to  maximum  orbital  velocity.  The  dashed  line  represents 
transition  to  turbulence  from  Figure  t.  The  solid  line  is  a  least  squares  fit 
of  an  expocential  curve  to  the  values  obtained  ai  the  onset  of  roiling 
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Fjg.  10.  Specwm  with  24  DOF  of  hoi  film  velocities  composed  from  *n 
ensemble  of  three  different  512-s  mos.  For  reference,  the  dashed  line  >s 
ensemble  avenge  (4£  DOF)  of  six  spear*  of  cross-shore  veJocroes  sho*  a 
in  Figure  3. 

A  power  spectrum  with  24  degrees  of  freedom  (DOF)  formed 
from  an  ensemble  of  three  hot  film  spectra  ts  show  n  in  Figure  10 
These  spectra  were  computed  from  data  runs  which  were  entirely 
unaffected  by  kelp  or  other  types  of  biofouling  and  were  fully 
dereciifted  pnor  to  transforming.  Inspection  shows  that  the 
spectral  rolloff  below  1  Hz  follows  the  Kolmogorov  -5/3  power 
law.  supporting  the  possibility  of  a  Turbulent  inertial  range  While 
the  turbulence  in  these  boundary  layers  is  intermittent,  this  is  not 
inconsistent  with  the  -5/3  spectrum  j5he  and  Orszcg.  1991).  In 
an  earlier  report  \lrman  ci  aJ..  1986].  it  was  reported  that  ibe  hot 
film  data  roll-off  bad  a  -2  slope,  but  w-e  now  believe  this  to  be  an 
artifact  resulting  from  recuficaiion.  The  bigber  frequencies  in 
Figure  10  are  polluted  by  what  is  believed  to  be  the  natural 
frequency  of  tbe  frame  of  the  surf  spider  (42  Hz)  and  by  noise 
from  the  60- cycle  power  tource. 

Using  the  results  cf  Kim  «  at  1 1971),  it  can  be  shown  (e  g.. 
Hinze.  1975)  that  for  steady  turbulent  flow,  the  turbulent  burst 
period  Th  scales  as 


where  5  is  the  boundary  layer  thickness.  Now  if  a  flow  over  a 
quasi-passive  tracer  is  experiencing  turbulent  bursting  such  that 
fluid  originating  at  tbe  tracer  interface  is  periodically  being  thrust 
into  the  free  stream,  one  might  expect  a  visual  pattern  of  puffs  of 
tracer  with  a  stream  wise  spac:ng  lh,  where 

lk  -  uTt  -  56  (9) 

Notice  how  similar  ibis  is  to  ibe  earlier  calculated  wavelength  for 
the  most  unstable  T-S  wave.  /t  -  76.  Using  the  boundary  layer 


thicknesses  of  3  cm  given  earlier  for  the  June  5.  1989.  experiment, 
equation  (9)  would  suggest  a  burst  spacing  of  *15  cm  Another 
indication  of  tbe  boundary  layer  thickness  would  be  ibe  tf.iLkriess 
of  the  granular- fluid  layer.  The  mean  thickness  of  ail  51  foiled 
layers  as  measured  from  the  at  rest  bottom  io  ibe  lop  of  :he  ru't 
surface  was  found  to  be  1  6  cm  with  a  standard  devotion  tf 
0.4  cn-  Using  this  as  the  boundary  layer  thickness  gnes  a  b-\r'! 
spacing  of  8  cm  While  these  estimates  of  the  turbulent  burst 
spacing  in  a  steady.  unidirectional  How  do  not  match  exactly  with 
ibe  measured  tuft  spacing  of  5-7  cm.  it  is  cerum  ihai  tbe  surface 
character  of  tbe  roiled  layer  suggests  a  process  dominated  t\ 
ejection  of  granular  laden  fluid  from  somew  here  tostde  the  reded 
layer. 

Tbe  mechanism  by  which  this  apparently  turbulent  flow 
develops  coherent  cross-stream  structures  as  evidenced  b>  the 
wales  in  tbe  granul.v-fluid  mixture  in  Figure  6  is  unknown 
However,  tbe  appearance  of  coherent  structures  in  turbulent  f)o«s 
is  not  a  new-  phenomenon  [e  g..  Mirujnt  end  Breland.  1974)  In 
visual  studies  of  a  turbulent  oscilla'ory  boundary  layer.  Hayasht 
and  Ohashi  [1982]  observed  coherent  large-scale  vertices  with 
cross-stream  axis  occurring  in  ibe  boundary  lay  er  just  prior  tc  the 
reversal  of  motion. 

Pluming 

Tbe  final  regime,  "pluming.*  also  occurs  under  rebeve’y  k  w 
amplitude,  long  period  w  aves  where  ibe  roiled  flow  ha*  developed 
a  periodic  structure  (Figure  4)  At  or  juM  after  ihc  passage  of  the 
wave  ctcsi  w ben  tbe  flow  begins  to  decelerate,  the  roiled  layer 
cxpJos/veh  hfis  off  ihc  bed.  injecting  aedi mom-laden  fluid  into  she 
interior  of  ihc  fluid  as  high  as  5-10  times  ihc  roiled  layer 
thickness  (Figure  II).  The  sediment  eject  a  is  then  carried  by  the 
free  stream  until  falling  out  at  the  cessation  of  motion,  leaving 
once  again  a  fiat  bed. 

The  mechanism  involved  in  pluming  is  difficult  to  identify,  not 
only  because  of  its  complexity  but  also  because  of  the  stringent 
conditions  apparently  required  for  its  occurrence  In  fact  no 
examples  of  pluming  were  observed  id  ihc  June  5  experiment,  so 
onset  criteria  bad  io  be  csnmaicd  from  ear  her  experiments 
Nonetheless,  as  pluming  w  as  generally  observed  on  the  same  days 
as  tbe  coherent  wale -like  structures  in  roiling,  it  is  likely  that  ihcy 
*re  related.  Most  probably,  pluming  represents  the  violent  breakup 
of  vortices  in  the  advent  of  an  adverse  pressure  gradient  associated 
with  the  passing  of  tbe  wave  crest.  The  longitudinal  vortices  with 
CTDss-sueaiD  axis  generated  during  the  reorganization  of  roiling 
appear  to  be  pressure  sensitive  and  experience  breakdown 
following  tbe  advent  of  deceleration. 

Relative  Occurrence 

Tbe  sequence,  streaking,  roiling,  and  pluming,  represents  an 
idealized  sequence  which  can  be  seeo  undcT  certain  conditions  but 
more  commonly  is  present  only  in  pan.  For  exampie,  undeT  very- 
low  wave  conditions  or  deeper  water,  the  only  poruon  which  may 
exist  would  be  streaking.  Under  more  general  conditions  (as 
exemplified  by  »he  June  5  experiment)  with  intermediate-period 
waves  or  incident  waves  with  a  broad  spectral  distribution,  the 
development  generally  only  includes  streaking  and  roiling  where 
tbe  bottom  is  flat  at  the  reversal  of  motioo.  This  "mild"  sequence 
is  schematized  by  case  A  in  Figure  12  In  relatively  "dean*  swell 
the  development  might  be  moderate  including  streaking  and  fully 
developed  roiling  resulting  in  a  bed  with  transition  npples  (case  B. 
Figure  12).  Under  relatively  low  amplitude,  long  period  swell 
where  wale-hke  structures  in  roiling  are  well  formed,  the  fully 
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developed  sequence,  streaking -roiling -pluming,  occu:?  re 
a  flat  bed  (case  C.  Figure  12i  A  final  case  under  pa 
intense  flows  character.?^  by  a  rapidly  growing  gran 
layer  is  shown  as  case  D  in  Figure  12  This  case  may  cr. 
portion  of  all  the  previous  cases  but  occurs  too  rap:, 
dissected  by  the  present  methods 


gradual  process  While  h*w  steepness  ripples  have  beer  observed 
under  :nc  trough  of  the  wave,  thus  arc  the  result  cf  waxes  with 
low  Shields  number  and  were  never  observed  tr  form  immediately 
foil'- w, rg  sheet  flow-  No  plum  it  g  has  been  observed  under  the 
wave  trough  in  die  field 

LABORATORY  OPMIRVaTIOsa 


Boundary  La yer  Development  under  the  Wave  Troyc.h  r.  ,,  ,  .  . 

Following  the  held  observations.  «»n  attempt  was  made  to 

The  observations  of  boundary  layer  development  under  the  wave  recreate  these  granular  fluid  flow  regimes  »c  the  laboratory  A  bed 
trough  show  a  sinking  contrast  to  that  under  the  crest.  The  initial  of  loose  sand  5  cm  deep  and  225  cm  long  was  placed  in  a  -a-:  5  rn 
stages  of  the  boundary  layer  development  under  the  trough  appear  long  by  200  cm  w  ide  wave  channel  The  surf  spider,  complete 
to  be  qualitatively  the  same  as  those  under  the  crest,  starting  v  :th  with  the  same  instrumentation  as  in  the  field  studies  minus  the 
a  laminar  boundary  layer  which  lead*  sr»tc  the  on -c;  cf  me  t*.  r.  video  vo^ueio.  was  mounted  m  me  center  of  the  bed  Shallow 


followed  by  streaking.  The  subsequent  regimes  of  roiling  and 
plumiDg  were  never  observed  to  develop  under  ibe  wave  trough 
Indeed,  development  of  the  granular-fluid  layer  under  the  trough 
appears  to  consist  instead  of  a  steady  thickening  of  a  granular-fluid 
sheet  without  the  abrupt  transitions  to  roiling  and  pluming  The 
surface  of  this  sheet  flow  exhibits  random  fluctuations  in 
cross-stream  elevation,  but  the  elevation  in  the  stream  wise 
direction  is  more  or  less  uniform.  In  contrast  to  tbe  roiied  layer, 
the  surface  of  this  sheer  suggests  the  predominance  of  turbulent 
mixing.  Tbe  thickness  of  the  granular-fluid  layer  for  similar  free 
stream  velocities  appears  to  be  smaller  under  the  trough  than  under 
tbe  crest  It  is  not  implied  that  transition  to  turbulence  docs  not 
occur  under  tbe  trough  On  the  contrary  ,  the  visual  impression  is 
that  the  granular- fluid  sheet  under  intense  conditions  is  most 
clearly  turbulent.  However,  there  is  no  clear  demarcation  for  the 
transition  to  turbulence,  and  it  appears  that  the  transition  is  3 


water  waves  of  5-s  periods  and  5 1 -cm  wave  heights  were 
generated  in  water  ujih  a  depth  c  l  132  cm  The  bed  was  filmed 
with  a  movje  camera  viewing  through  a  Plexiglas  window  whije 
synchronized  measurements  were  recorded  as  in  die  field 
experiments  The  roam  observation  from  these  tests  was  that  no 
roiling  or  pluming  occurred  and  that  any  as  vm  me  try  between  crest 
and  trough  was  not  readily  apparent  or  could  be  attributed  to 
velocity  and  acceleration  asymmetries  inherent  in  the  shallow 
water  wave  form  This  is  in  contrast  to  our  field  experiments 
where  waves  with  a  similar  Reynolds  number  \  He  -  ?  5  x  10  ) 
clearly  developed  into  roiling  (Figuic  ui  and  exhibited  the 
associated  crest  trough  asymmetry 

DlSCt'SSION 

Other  investigators  je  g  .  Hants  an,i  f-Juntlcx.  i'-br  j  have 
observed  crest  trough  asymmetries  in  the  sediment  nsptir.se  and 
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Fjg.  12.  Schematic  drawing  of  four  identifiable  sequences  in  the 
development  of  the  fluid -granular  boundary  layer  under  the  crest  of  near- 
breahng  waves,  shown  ia  the  order  of  increasing  intensity  of  motion  at 
the  interface.  Details  of  case  C  are  shown  in  Figure  4. 


have  c'lributed  ibis  behavior  to  asymmetries  in  the  free  stream 
velocity  and/or  acceleration.  However.  Figure  13  indicates  that 
such  an  explanation  cannot  fully  explain  (he  asymmetry  observed 
in  our  experiments.  A  histogram  of  (be  maximum  onshore  (crest) 
and  offshore  (trough)  velocities  observed  during  the  June  5 
experiment  is  shown  in  Figure  1 3a  together  with  the  observed 
range  of  free  stream  velocities  at  the  onset  of  roiling  Clearly,  the 
offshore  flows  are  not  velocity  limited,  yet  roiling  did  not  occur. 
In  Figure  13b  we  have  plotted  the  maximum  onshore  Eulerian 
accelerations  as  well  as  the  observed  range  of  accelerations  at  the 
onset  of  roiling.  For  the  offshore  motion,  we  have  plotted  the 
observed  trough  acceleration  maxima  for  only  the  wave  half -cycles 
which  had  free  stream  velocities  within  the  observed  raoge  for  the 
onset  of  roiling.  Again  it  is  clear  that  the  trough  flow,  when 
compared  to  the  crest,  is  not  acceleration  limited. 

Many  theoretical  investigations  (e  g.,  Longuet-Higgins,  1953; 
Jacobs,  1984]  predict  a  steady  mass  transport  in  the  boundary 
layer  beneath  progressive  gravity  waves.  While  observations  of 
this  phenomenon  have  been  made  in  the  laboratory,  there  was  little 
mean  boundary  layer  current  observed  in  our  near-breaker 
experiments.  The  three  hot  Film  runs  used  to  compute  the 
ensemble  spectrum  which  appears  in  Figure  10  bad  probe  heights 
of  13  cm.  4.2  cm.  and  1.8  cm  above  the  at -rest  bed  pnor  to  the 
start  of  each  run  The  Field  calibrations  for  these  three  runs  had 


standard  deviations  with  respect  to  the  em  current  meter  of  0  9%. 
i.2%,  and  0.9%  respectively,  where  the  values  are  expressed  as  3 
percentage  of  the  measured  values  The  mean  current  ever  each 
of  tbc  three  data  runs  was  0  5.  2.5.  and  1.9  cm/s  respectively  It 
is  clear  from  Figure  13  that  meaa  currents  of  this  magnuude  are 
not  sufficient  to  explain  tbc  asymmetry  observed  in  our 
experiments. 

Numerous  theoretical  investigation*  of  progressive  graw rv 
w  ave:  propagating  over  a  permeable  bed  [Reid  and  Kajiura.  195'  . 
Packwood  and  Peregrine.  1980)  indicate  that  an  osculauog  flow 
is  induced  in  the  bed  and  results  io  flow  normal  to  the  bed  surface 
under  the  crest  and  trough  of  the  wave.  This  flow  is  the  result  of 
the  spatially  varying  pressured  Field  imposed  on  the  bed  over  an 
entire  wavelength.  This  ventilation  flow  is  out  of  the  bed  and  into 
the  fluid  (injection)  under  the  wave  trough  and  into  the  bed 
(suction)  under  tbe  crest  These  studies  also  indicate  that  with 
typical  natural  beds,  the  induced  velocities  are  several  orders  of 
magnitude  smaller  than  the  horizontal,  free  sueam.  wave-orbital 
velocities,  and  so  any  effect  on  tbe  free  stream  velocities  is  small 
On  tbe  other  band,  tbe  effect  of  small  normal  flows  on  the 
oscillatory  boundary  layer  is  not  necessarily  negligible. 

It  has  long  been  known  (see  for  example  Tewfik  (19631.  Mtckley 
and  Davis  [19571.  Schlichiing  (1979ji  that  except  for  transitional 
processes,  boundary  sucuoa  with  steady  flow  leads  to  enhanced 
bed  stress  and  injection  leads  to  reduced  bed  stress  even  for  the 
small  ratios  of  boundary  normal  to  free  stream  velocity  u>  be 
expected  in  in  oscillatory  boundary  layer  over  a  sand  bed  It  is 
also  known  that  suction  has  tbe  effect  of  stabilizing  the  flow, 
reducing  the  three-dimensionality  of  tbe  turbuleot  boundarv  layer 
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MAX  ACCELERATION  (cw/w1) 

Fig.  13.  Histograms  of  (a)  orbital  velocity  acd  (b)  acceleration  maxima 
associated  with  onshore- offshore  orbital  motion  for  the  Field  deployment 
of  Juoe  5.  1989.  Solid  brackets  indicate  ranges  of  observed  roiling  onset 
during  onshore  (crest)  mouon.  Dashed  brackets  indicate  same  ranges  for 
offshore  (trough)  motion. 
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and  increasing  the  relative  coombutioo  of  large-scale  coherent 
structures  [Fulachier  el  at.  1982}.  Tenfik  (19631  reports  that 
injection  decreases  the  rate  of  growth  of  the  turbulent  boundary 
layer.  Injection  also  moves  the  shear  maxima  away  from  the  bed 
which,  combined  with  a  “law  of  the  wall'  type  approximation, 
suggests  that  under  the  trough  sediment  suspending  turbulent 
kinetic  energy  (TKJE)  is  removed  from  the  bed  w  hich  is  the  source 
of  sediment  The  reverse  is  true  of  suction. 

It  is  suggested  from  the  above  that  an  oscillatory  boundary  layer 
over  a  permeable  satrd  bed.  alternating])'  experiencing  suctioo  and 
injection  under  the  crest  and  trough  of  the  wave  respectively, 
would  exhibit  the  following  characteristics.  The  turbulent 
boundary  layer  under  the  trough  would  be  expected  to  be  slow- 
developing  aod  fully  homogeneous  with  little  structure.  Any 
internal  dynamics  within  the  boundary  layer  under  the  trough 
would  be  obscured  by  strong  turbulent  mixing  which  would  occur 
further  from  the  bed.  around  the  displaced  TKE  maxima.  The 
sediment-mo bilaing  properties  of  the  flow  would  be  dummsbcd 
(decreased  bed  stress  with  TK£  removed  from  the  bed)  which 
would  be  reflected  in  thinner,  less  dense  graDular-fiuid  layers. 
Flow  under  the  crest  would  be  characterized  by  a  more  rapid  aod 
therefore  distinct  turbulent  boundary  layer  development.  The 
boundary  layer  would  be  expected  to  be  less  homogeneous, 
exhibiting  more  structure  than  under  the  trough.  Sediment 
mobilization  would  be  enhanced.  Such  a  phenomenon  would 
clearly  have  strong  repercussions  for  any  type  of  modeling  wfcucb 
treats  boundary  layer  development  as  a  function  of  magnitude 
only.  For  example,  the  concept  of  friction  factors  m  this  type  of 
5)iuai20D  would  require  a  separate  set  of  friction  factors  for  flow 
UDder  the  crest  and  trough  respectively. 

It  should  be  noted  that  turbulent  transition  under  the  crest  which 
occurs  at  a  phase  greater  than  90*  would  be  experiencing  the 
destabilizing  affect  of  an  adverse  pressure  gradient  and  would  be 
expected  to  be  qualitatively  different  than  transition  which  occurs 
earlier.  eveD  though  suction  could  still  be  playing  a  roll  m 
boundary  layer  character! sties. 

For  permeable  beds  where  the  thickness  of  the  bed  is  small 
relative  to  the  wavelength.  Pachvood  and  Peregrine  [1980]  shew 
that  the  ratio  of  vertical  velocity  to  free  stream  velocity.  H 
under  the  crest  of  the  wave  is  linearly  proportional  to  the  thickness 
of  the  permeable  bed.  Thus  the  crest-trough  asymmetry  observed 
in  tbe  field  would  not  be  expected  to  be  seen  in  a  laboratory 
setting  which  does  not  have  a  permeable  bed  which  is  either 
relatively  deep  and/or  Iodj  enough  to  simultaneously  “feel"  an 
entire  wavelength. 

Aii  >ui  observations  support  this  theory  aod  a  laboratory 
experiment  utilizing  a  venlilaied  bed  is  under  way  in  order  to 
investigate  this  phenomenon  more  closely. 

Conclusions 

In  situ  field  observations  have  been  made  of  the  development 
of  tbe  wave-induced  fluid-granular  boundary  layer  over  penneab’e 
beds  of  loose  saDd.  The  visual  record  accompanied  by 
synchronized  measurement  of  physical  parameters  leads  to  the 
following  conclusions: 

1.  The  full  fluid-granular  boundary  layer  development  under 
tbe  crest  of  tbe  wave  can  be  divided  into  an  identifiable  three-pan 
sequence  composed  of  streaking,  roiling,  and  pluming. 

2.  Tbe  streaking-roiling  transition  is  tbe  manifestation  of  a 
sudden  laminar -turbulent  transition. 

3.  A  crest-trougb  asymmetry  m  tbe  fluid-granular  boundary 
layer  development  has  been  identified  This  asymmetry  is 


independent  of  asymmetries  id  the  near  bed  orbsUl  velocity  and 
acceleration  and  is  thought  to  be  related  ic  effects  of  veclilauon 
m  the  permeable  bed. 
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m  the  terminology  of  sequences  io  the  development  of  the  fluid  granular 
boundary  layer. 
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CHAPTER  II 


VENTILATED  OSCILLATORY  BOUNDARY  LAYERS 
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1.  INTRODUCTION 

Laboratory  studies  of  the  boundary  layer  arising  from  oscillatory  flow  have  always 
treated  the  problem  as  a  cyclical  process  composed  of  two  purely  symmetrical  half  cycles  [eg 
Jonsson,  1963;  Sleath  1987],  Asymmetries  occurring  in  the  laboratory  [Rick  ct  al,  1981]  or  the 
field  [Hanes  &  Huntley,  1986]  have  been  explained  by  asymmetries  in  the  fluid  velocity  or 
acceleration  inherent  in  shoaling  waves.  In  a  paper  reporting  on  field  observations  of  w'ave  driven 
oscillatory  boundary  layers  over  sand  beds,  Conley  and  Inman  [1992]  reported  the  observation  of 
an  asymmetry  in  the  development  of  the  boundary  layer  which  could  not  be  explained  by  velocity 
or  acceleration  asymmetries  in  the  overlying  fluid.  As  a  possible  mechanism  for  this  asymmetry, 
they  proposed  a  type  of  transpired  boundary  layer  which  is  here  called  the  ventilated  oscillatory 
boundary  layer.  This  work  reports  on  laboratory  experiments  designed  to  investigate  the  effect  of 
boundary  ventilation  on  the  velocities,  bed  stress  and  turbulent  flow  properties  of  the  turbulent 
oscillatory  boundary  layer. 

Transpired  boundary  layers  have  practical  significance  in  many  applications  and  have 
been  the  subject  of  studies  for  well  over  forty  years  [eg  Libby  et  al,  1952;  Antonia  et  al,  1990]. 
These  arc  boundary  layers  which  arise  from  a  fluid  flowing  with  velocity  u  parallel  to  a  solid 
permeable  boundary  through  which  a  secondary  fluid  is  flowing  with  velocity  w.  The  secondary 
flow  through  the  boundary  is  generally  referred  to  as  transpiration.  If  the  sense  of  transpiration  is 
into  the  overlying  fluid  (w  >  0),  it  is  called  injection  or  blowing.  If  the  sense  of  transpiration  is 
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from  the  fluid  into  the  surface  (vv  <  0)  it  is  termed  suction.  In  general  the  two  fluids  are  identical 
and  tiic  transpired  flow  is  smaller  in  magnitude  than  boundary  parallel  llow  so  that  the  transpira¬ 
tion  parameter  V=w/u  has  a  magnitude  less  than  1. 

Ventilated  oscillatory  boundary  layer  is  a  term  applied  here  to  the  boundary  layer  arising 
from  an  oscillating  boundary  parallel  flow  subject  to  oscillating  transpiration  of  the  same  period 
and  shape.  In  this  work,  the  dominant  flow  is  of  the  form  u(t)-  umsin  (or)  and  subject  to  tran¬ 
spiration  of  the  form  w(t)  -  wmsin(at+$).  Here  o  =  2n/T  is  the  radian  frequency,  T  is  the 
period  and  <{>  is  an  arbitrary  constant  phase.  Analogous  to  the  transpiration  parameter  V',  a  flow 
ventilation  parameter  V  shall  be  defined  as  V  =  wm/um.  The  sign  of  V  indicates  whether  injec¬ 
tion  (V >0)  or  suction  (V<0)  occurs  concurrently  with  positive  flow.  Additionally,  an  instantane¬ 
ous  ventilation  parameter  shall  be  defined  as  V"  =  w(t)/u(t).  Through  most  of  this  work,  the 
phase  0  will  be  0  in  which  case  V'  becomes  a  constant,  identically  equal  to  V. 

In  general  boundary  transpiration  affects  the  boundary  layer  velocity  profile.  Suction 
tends  to  pull  streamlines  down  closer  to  the  bed  shifting  the  velocity  profile  closer  to  the  wall. 
This  results  in  higher  shear  near  the  bed  and  therefore  higher  shear  stress  at  the  bed.  This  is 
demonstrated  quite  dearly  in  a  comparison  of  the  asymptotic  suction  profile,  an  exact  solution  to 
the  Navier  Stokes  equations  for  flow  over  a  flat  plate  with  constant  suction,  and  the  Blasius 
profile  [Schlichting,  1979).  In  contrast,  injection  results  in  a  spreading  of  the  streamlines  near  the 
boundary  which  reduces  the  near  bed  shear  and  the  resultant  bed  stress. 

In  a  study  of  turbulent  boundary  layers  in  steady  flow  over  flat  plates  with  injection, 
Mickley  &  Davis  [1957]  show  that  the  friction  factor  Cf  is  a  strong  function  of  V  and  Re*  where 
Re*  is  the  Reynolds  number  based  on  distance  along  the  plate.  However  their  results  can  be 
shown  to  indicate  that  the  ratio  Cf/Cfo  where  Cfo  is  the  friction  factor  with  out  injection,  is  a 
strong  function  of  V  but  exhibits  almost  no  Re*  dependency.  The  results  of  Simpson  et  al  [  1969] 
show  the  same  to  be  true  for  steady  suction  on  a  fiat  plate. 


The  asymptotic  suction  profile  shows  reduced  curvature  in  the  boundary  layer  velocity 
profile  relative  to  the  Blasius  profile  [Schlichting,  1979).  This  suggests  that  suction  tends  to  sta¬ 
bilize  the  flow.  Analytical  solutions  for  boundary  layers  with  injection  show  an  inflection  point 
in  the  velocity  profile  suggesting  that  injection  tends  to  destabilize  the  flow.  Experiments  have 
shown  these  findings  to  be  true  not  only  for  laminar  flow  but  turbulent  flow  as  well.  Tewfik 
[1963]  demonstrated  that  turbulent  boundary  layers  with  injection  were  thicker  than  would  bo  the 
case  without  injection.  It  was  also  shown  that  while  injection  reduces  stress  at  the  bed  and 
through  out  the  inner  tenth  of  the  boundary  layer,  it  increased  it  elsewhere.  Antonia  et  al  [1988] 
reconfirmed  that  suction  results  in  a  significant  reduction  in  turbulent  velocity  fluctuations  u'  and 
w'.  While  investigating  the  effect  of  suction  on  organized  motion  in  turbulent  boundary  layers, 
they  observed  that  low  speed  streaks  in  the  near  wall  region  were  more  persistent  and  that  the  fre¬ 
quency  of  dye  injections  into  the  outer  layer  was  reduced.  This  is  consistent  with  the  results  of 
Fulachier  [1982]  where  the  rate  of  turbulent  production  was  found  to  decrease  and  the  relative 
contribution  of  large  scale  structures  was  found  to  increase. 

Kays  [1972]  shows  that  transpired  turbulent  boundary  layers  seem  to  respond  to  local 
conditions.  That  is  a  boundary  layer  experiencing  a  sudden  change  in  injection  will,  following  a 
short  period  of  adjustment,  behave  as  though  it  had  always  been  subject  to  that  level  of  injection. 
This  appears  to  be  less  true  for  a  reduction  in  injection  than  for  an  increase.  It  is  also  shown  that 
transpired  boundary  layers  are  sensitive  to  acceleration,  the  sensitivity  increasing  with  V7.  No  stu¬ 
dies  have  been  reported  on  the  effects  of  transpiration  in  oscillatory  flow. 

2.  EXPERIMENTAL  SETUP 

The  experiments  have  been  performed  using  the  oscillatory  flow  tunnel  (OFT)  at  the 
Scripps  Institution  of  Oceanography  hydraulics  laboratory  (Figure  1).  The  OFT  is  composed  of 
an  acrylic  working  section  7.0  m  long  which  terminates  in  painted  steel  cylindrical  risers  at  each 
end.  One  riser  acts  as  a  reservoir  and  is  open  to  the  atmosphere,  the  other  riser  is  scaled  by  a  62 
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cm  diameter  acrylic  piston  which  drives  the  oscillating  How.  The  piston  is  controlled  by  a 
hydraulic  ram  which  can  be  programmed  to  respond  to  any  electrical  input.  Further  details  of  the 
facility  arc  described  by  King  et  al  [1984],  For  these  experiments  a  false  floor  was  placed  in  the 
working  section  giving  a  tunnel  cross  section  39  cm  wide  and  28  cm  deep.  The  test  section  was 
created  by  placing  a  rigid,  permeable  floor  in  the  center  4.25  m  of  the  working  section.  This  floor 
was  composed  of  fourteen  30  x  30  x  2.5  cm  bricks  of  fused  Alundum.  These  bricks  were  com¬ 
posed  of  Alundum  grains  with  a  mean  diameter  of  250  microns  which  were  fused  under  high 
temperature.  Prior  to  fircing,  the  bricks  were  pressed  at  high  pressure.  As  a  result,  the  flat  sur¬ 
faces  of  the  rather  angular  grains  aligned  to  form  a  smooth  surface  on  the  bricks.  Roughness  ele¬ 
ments  arc  caused  by  pores  in  the  surface,  not  by  protuberant  grains.  The  mean  pore  size  of  the 
internal  pores  is  stated  to  be  180  microns.  A  secondary  piston  which  was  plumbed  into  the  the 
cavity  below  the  test  section  was  used  to  drive  the  boundary  layer  ventilation. 

Tests  were  performed  using  a  period  of  oscillation,  T  =  7  s  and  an  orbital  displacement 
d o  of  207  cm.  This  results  in  a  flow  with  oscillatory  Reynolds  number.  Re,  of  l.OxlO6  where 
Re  is  defined  by 


Re  = 


“mdo 

2v 


(1) 


where  um  is  the  maximum  mass  averaged  velocity.  The  results  of  Jensen  ct  al  [  1989]  have  shown 
that  oscillatory  boundary  layers  of  this  Re  will  be  fully  turbulent  for  phases  greater  than  45°  and 
transitional  for  phases  greater  than  15°.  For  the  unventilated  case,  the  boundary  layer  thickness 
8q  is  1.8  cm  where  5o  has  been  defined  as  the  first  level  at  which  the  vertical  derivative  of  velo¬ 
city  vanishes  at  a  phase  of  90°.  The  maximum  friction  velocity  u*m  =  ^Tm  / p  is  4.1  cm/s.  The 
N'ikuradse  equivalent  roughness  pirameter  ks  as  determined  from  boundary  layer  velocity 
profiles  was  found  to  be  0.003  cm  giving  a  roughness  Reynolds  number  k^-u*mkslv  of  1.3.  The 
ventilation  parameter  V  was  zero  in  the  control  tests  and  ranged  in  magnitude  from  2.5x1 0-4  to 
4.x  10  in  other  tests.  The  amplitude  of  w  was  calculated  by  assuming  uniform  velocity  over 


the  test  bed  and  distributing  the  volumetric  displacement  of  the  secondary  piston  over  the  enure 
test  bed.  Calculations  indicate  that  uniform  velocity  over  the  test  bed  is  a  good  assumption  gi  ven 
the  manufacturers  supplied  permeability  of  192  Darcys  for  the  Alundum  bricks.  This  assumption 
was  tested  through  ^ye  studies  in  which  no  systematic  variations  in  the  How  through  the  bricks 
could  be  detected. 

Three  main  flow  parameters  are  presented  in  this  work.  These  include  the  bed  stress 
X(t),  fluid  velocity  u(z,t)  and  the  vertical  component  of  fluctuating  velocity  The  bed 

shear  stress  was  measured  using  a  TSI  model  1237W  flush  mounted  hot  film  sensor  as  shown  in 
Figure  2.  The  sensor  was  locked  into  an  acrylic  probe  holder  9.5  mm  in  diameter  which  fit 
through  an  O-ring  seal  into  a  cylindrical  acrylic  sleeve  of  1.8  mm  wall  thickness.  The  sleeve  was 
permanently  epoxied  into  the  test  bed  with  the  top  of  the  sleeve  flush  with  the  tops  of  the  rough¬ 
ness  elements  forming  the  fixed  surface  of  the  Alundum  bed.  When  inserted,  the  tops  of  the 
probe  and  holder  were  level  with  the  top  of  the  sleeve.  This  arrangement  allowed  the  shear  stress 
sensor  to  be  calibrated  in  a  separate  facility  and  maintain  the  same  relationship  between  probe 
and  holder  when  transferred  to  the  test  facility. 

The  longitudinal  component  of  velocity  was  measured  using  a  TSI  model  1210-20W 
cylindrical  hot  film  sensor.  The  probe  was  mounted  on  a  sting  projecting  through  the  lid  of  the 
OFT  with  the  axis  of  the  sensing  element  horizontal  and  transverse  to  the  direction  of  flow.  This 
probe  was  calibrated  by  placing  it  in  the  tunnel  center  and  calibrating  agaiast  the  tunnel  mass 
averaged  velocity  using  a  technique  adopted  from  Flick  &  George  [1990]  While  the  hot  film 
measures  fluid  speed,  the  time  series  was  dercetified  using  the  mass  averaged  fluid  velocity  zero 
crossings  as  a  reference.  Turbulent  velocity  fluctuations  were  measured  using  a  TSI  model 
1287W  split  film  hot  film.  This  sensor  has  a  cylindrical  sensing  element  composed  of  two  halves 
which  arc  electrically  insulated.  This  design  permits  decomposition  of  the  measurements  into 
two  components  of  velocity.  Briefly,  the  magnitude  of  the  cooling  for  both  halves  is  dependent  on 
the  magnitude  of  the  velocity  vector  and  the  ratio  of  the  cooling  is  dependent  on  angle  the  vector 
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makes  with  the  plain  of  the  splits.  Operation  of  the  sensor  is  discussed  in  [Blinco  &  Sandbom, 
1975].  The  directional  sensitivity  was  calibrated  once  before  and  after  all  tests  while  the  velocity 
magnitude  response  was  calibrated  for  every  test  as  described  for  the  cylindrical  hot  films  above. 


Depending  on  die  test,  data  runs  were  of  500  or  700  seconds  in  duration.  The  time 
scries  collected  over  each  run  were  ensemble  averaged  in  order  to  determine  the  phase  averaged 
quantities.  In  oscillatory  flows,  the  turbulent  decomposition  of  a  flow  property  %(t)  has  three 
components,  the  time  mean  <^>,  the  easemble  averaged  mean  ^(0)  and  the  ensemble  averaged 

r.m.s.  mean  of  the  fluctuating  quantity  (0).  In  a  purely  oscillatory  flow,  the  traditionally- 
defined  time  mean  is  identically  zero.  For  a  time  series  collected  over  N  cycles  of  period  T,  the 
ensemble  average  mean  is  calculated  as 


1  H  f  "1 

£(or)  =  —  «  +  (t-])7~ 

/vi=  t  1  J 


(2) 


2n 

where  o  is  the  radian  frequency  ~jr  and  the  definition  is  valid  for  t~ 0  to  t~T.  Similarly  the 

ensemble  averaged  r.m.s.  mean  of  the  fluctuating  quantity  is  defined  by 

i  n  r 

*  * 


k’2(ot)  = 


N- 1 


ct  +  ( i-\)T 


-^(cr) 


(3) 


This  decomposition  is  schematicized  in  Figure  3.  Slcath  1 1987)  and  others  have  shown  that  there 
was  no  significant  im  ovement  in  the  consistency  of  the  estimate  of  ensemb’e  averaged  mean 
quantities  for  averages  of  greater  tha^  50  cyJes.  The  500  s  tests  in  this  work  represent  71  cycles 
and  the  700  « tests  represent  100  cycles.  When  discussing  the  phase  of  the  oscillatory'  flow  0,  the 
convention  will  be  used  that  9=0  corresponds  to  the  negative  to  positive  velocity  zcio  crossing. 
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3.  EFFECTS  OF  VENTILATION 

The  ensemble  averaged  centerline  velocity  for  the  ease  of  no  flow  through  the  bed  is 
shown  in  Figure  4a.  As  can  be  seen  from  the  plot,  the  velocity  is  relatively  symmetric  throughout 
the  entire  cycle.  Each  half  cycle  appears  qualitatively  similar  tc  the  other  ha'f  cycle  and  the 
accelerating  portions  of  each  half  cycle  appear  to  be  mirror  images  of  the  decelerating  ponion. 
Figure  4b.  is  a  plot  of  ensemble  averaged  centerline  acceleration.  This  plot  reveals  asymmetries 
in  the  flow  which  are  not  obvious  in  the  velocity.  In  particular  it  is  seen  that  acceleration  tends 
towards  a  monotonical  decrease  through  out  the  first  half  cycle  as  vould  be  expected  in  an  oscil¬ 
latory  flow.  This  is  not  true  for  the  second  half  cycle  where  periods  of  decreasing  acceleration  are 
clearly  present  in  what  would  be  expccied  to  be  a  time  of  monotonically  increasing  acceleration. 

While  these  fluctuations  don’t  appear  to  have  a  big  affect  on  fluid  velocity,  it  might  be 
anticipated  that  these  fluctuations  in  pressure  gradient  would  have  a  larger  affect  in  the  boundary 
layer  and  on  the  bed  stress.  This  seems  to  be  bom  out  by  the  results  presented  in  I  ,gure  4c.  The 
solid  line  represents  the  ensemble  averaged  bed  stress  as  derived  from  the  shear  stress  sensor. 
The  bed  stress  in  the  second  half  cycle  is  much  peakier  than  would  be  expected  in  an  oscillatory 
flow  of  this  Reynolds  number  Re  =  106  while  the  first  half  cycle  is  more  typical  [see  Jensen  et  al. 
1989], 

Jensen  el  al  [1989|  have  shown  that  the  boundary  layer  velocity  profiles  in  high  Rey¬ 
nolds  number  flows  do  exhibit  a  logarithmic-layer  from  which  the  Nikuradse  equivalent  rough¬ 
ness  iength  ks  and  the  bed  stress  can  be  determined.  This  was  found  to  be  especially  true  around 
the  occurrence  of  peak  free  stream  velocity.  In  order  to  test  this  in  the  present  experiments,  ks 
was  taken  as  the  average  value  derived  from  the  best  fits  to  a  logarithmic  layer  in  the  phase  range 
4S°<9>!  35°  This  value  of  ks  was  then  used  to  calculate  the  bed  stress  through  out  the  wave 
cycle.  These  values  are  plotted  as  heavy  dots  in  Figure  4c  and  as  can  be  seen,  they  give  a  reason¬ 
able  estimate  of  bed  stress  in  the  first  half  cycle.  However  when  this  same  procedure  is 


performed  using  the  roughness  derived  from  proliles  in  the  range  225°<0>3 15",  n  results  m  a 
underestimation  of  bed  stress  as  can  be  seen  by  the  open  dots  in  Figure  4c.  tn  fact  it  appears  that 
bed  stress  estimates  based  on  ks  derived  from  the  first  half  cycle  are  better  than  those  based  solely 
on  second  half  cycle  parameters.  The  logarithmic-layer  is  a  consequence  of  the  presence  of  a 
constant  stress  layer  in  the  boundary  layer.  It  is  not  surprising  that  the  presence  of  a  fluctuating 
pressure  gradient  would  disturb  this  constant  stress  layer  as  evidently  occurs  in  the  second  half 
cycle.  As  a  consequence  of  this  anomalous  behavior,  the  following  convention  is  adopted 
throughout  the  paper.  When  discussing  bed  stress,  results  from  only  the  first  half  cycle  shall  be 
used  and  the  sense  of  ventilation  will  be  alternated  to  provide  both  suction  and  injection  during 
that  half  cycle.  The  full  cycle  shall  only  be  utilized  when  presenting  velocity  and  turbulence 
profiles.  While  in  general  the  bed  stress  results  from  the  second  wave  half  cycle  differ  quantita¬ 
tively  from  the  first,  the  conclusions  arrived  at  in  this  work  are  qualitatively  supported  by  the 
results  from  both  half  cycles. 

Boundary  Layer  Velocity 

As  discussed  earlier,  the  most  direct  way  that  suction  is  expected  to  affect  the  boundary 
layer  is  to  draw  the  velocity  profile  down  close  to  the  bed.  In  the  boundary  layer  this  results  in 
higher  mean  velocities  near  the  bed  and  a  greater  shear  at  the  bed.  The  affect  of  injection  is 
expected  to  be  just  the  reverse,  namely  streamlines  pushed  away  from  the  bed  leading  to  lower 
mean  boundary  layer  velocities  and  a  reduced  bed  shear.  This  pattern  is  exactly  what  is  seen  in 
Figure  5.  This  figure  shows  a  series  of  boundary  layer  profiles  at  various  phases  throughout  the 
entire  cycle.  At  each  phase  there  are  two  profiles.  The  solid  line  represents  the  velocity  profile 
for  the  case  of  no  ventilation  and  the  dashed  line  represents  the  velocity  profiles  for  V  =  SxlCT4. 
It  is  immediately  obvious  that  even  for  such  a  small  value  of  the  ventilation  parameter,  the  mean 
velocities  in  the  boundary  layer  are  affected.  In  particular  it  is  seen  that  during  the  occurrence  of 
suction  (()°<0<18()H)  the  mean  velocities  throughout  the  boundary  layer  and  well  above  it  arc 
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uniformly  greater  than  the  unventilated  velocities.  This  fact  necessitates  that  the  shear  at  the  bed 
is  greater  which  is  in  fact  clear  in  the  figure.  Once  again,  just  the  reverse  is  true  for  injection 
(180°<9<360°).  It  is  even  more  readily  evident  that  the  mean  velocities  throughout  the  boun¬ 
dary  arc  lower  with  injection  than  the  unventilated  case  and  therefor  the  boundary  shear  is  less 
than  the  unventilated  case. 

Considering  the  symmetric  nature  of  unventilated  oscillatory  flow,  it  is  clear  that  the 
above  observations  require  an  asymmetry  in  the  mean  boundary  layer  velocities  for  a  ventilated 
oscillatory  boundary  layer.  In  fact  this  asymmetry  must  be  realized  as  a  net  current  in  the  boun¬ 
dary  layer  above  the  permeable  bed.  This  current  will  be  positively  directed  for  negative  V  and 
vise  versa. 


Turbulence 


A  contour  plot  of  the  development  of  vvv  over  phase  and  elevation  is  shown  in  Figure 
6.  This  is  for  the  case  of  no  ventilation  and  the  results  are  qualitatively  similar  to  those  reported 

by  others  [see  Sleath  1987J.  Sleath  [  1987]  reported  that  and  behave  qualitatively  the 


same,  so  for  the  following  discussion,  it  shall  be  assumed  that  the  behavior  of  vw'2  is  indica¬ 
tive  of  turbulent  intensities  in  general.  In  general  the  maximum  in  turbulence  intensity  is  seen  to 
occur  near  the  bed  just  following  the  start  of  deceleration  (90°  and  270°).  From  this  maxima  the 
turbulence  slowly  diffuses  upward  over  phase  giving  the  impression  of  an  inclined  plume.  For 
unventilated  oscillatory  flow,  this  pattern  is  relatively  symmetrical  with  turbulence  maxima 
occurring  at  about  the  same  level  in  each  half  cycle  and  the  siz.e,  shape  and  inclination  of  the  tur¬ 
bulent  plumes  are  the  same  in  each  case.  It  is  also  seen  that  the  turbulence  levels  associated  with 
each  half  cycle  return  to  the  background  level  prior  to  the  arrival  of  the  next  plume. 

The  effects  of  ventilation  on  the  turbulent  properties  of  the  flow  are  shown  in  Figure  7. 
This  is  a  contour  plot  of  the  development  of  over  phase  and  elevation  for  the  case  of  V  = 
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0.001 .  The  figure  shows  three  major  effects  of  ventilation.  The  first  effect  is  that  near  the  bed,  the 
turbulence  levels  arc  higher  with  suction  (0°<9<180°)  than  injection.  This  is  observed  in  the 

figure  where  the  contours  of  Vw'2  are  one  level  higher  for  the  suction  half  cycle  than  for  the 
injection  half  cycle.  Second,  with  suction  the  turbulence  maximum  is  drawn  closer  to  the  bed. 
Finally,  due  to  the  flow  stabilizing  tendency  of  suction  and  the  destabilizing  influence  of  injec¬ 
tion,  the  turbulent  plumes  have  become  totally  distorted.  In  particular,  the  turbulent  plume  ori¬ 
ginating  in  the  half  cycle  with  suction  is  very  much  diminished  and  is  contained  in  a  thin  layer 
(=5mm)  drawn  near  the  bed.  The  turbulent  plume  from  the  injection  half  cycle  is  much  enhanced, 
extending  higher  into  the  flow  and,  consequentially,  much  later  in  the  flow.  This  tendency  is  so 
pronounced  that  it  would  appear  that  the  majority  of  the  turbulence  that  arises  from  the  injection 
half  cycle  occurs  during  the  subsequent  half  cycle.  The  net  affect  of  this  is  that  a  "snapshot"  of 
turbulence  levels  at  constant  phase  would  show  the  highest  turbulence  levels  to  occur  during  the 
half  cycle  with  suction! 

Using  a  logarithmic  elevation  scale,  the  contour  plot  of  the  time  history  of  &  (Fig¬ 
ure  8)  shows  that  the  previously  described  qualitative  picture  is  true  for  V  as  low  as  -0.0005,  the 
lowest  ventilation  parameter  for  which  data  of  this  type  was  collected.  In  short,  the  turbulence 
levels  near  the  bed  are  enhanced  for  suction,  the  turbulent  plume  associated  with  suction  is 
confined  close  to  the  bed  and  the  turbulent  plume  associated  with  injection  is  enhanced  and  per¬ 
sists  in  time  throughout  most  of  the  subsequent  half  cycle. 

Bed  Stress 

The  effect  of  boundary  ventilation  on  bed  stress  is  summarized  in  Figure  9.  Plotted  in 
this  figure  is  the  ensemble  averaged  bed  stress  from  500  second  runs  for  various  ratios  of  the  ven¬ 
tilation  parameter  V.  The  heavy  solid  line  is  the  bed  stress  for  the  case  of  no  ventilation  and 
serves  as  the  baseline  for  the  other  tests.  The  peak  value  from  this  curve  has  been  used  to 
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normalize  all  the  data.  As  can  be  seen  in  the  figure,  ventilation  clearly  affects  the  bed  stress  in 
oscillatory  flow  with  suction  (V <0)  leading  to  increased  bed  stress  and  injection  (^>0)  causing 
reduced  bed  stress.  This  figure  shows  that  there  is  a  significant  stress  reduction  or  enhancement 
for  ventilation  parameters  with  magnitudes  as  low  as  2.5X10-4. 

The  data  from  Figure  9  has  been  used  to  construct  Figure  10.  For  each  value  of  V,  the 
bed  stress  was  integrated  over  all  positive  values  and  normalized  by  the  integrated  bed  stress  for 
no  ventilation.  These  points  are  plotted  as  a  function  of  ventilation  parameter  and  are  represented 
by  the  open  points  in  Figure  10.  The  solid  line  is  a  composite  of  two  functions  which  seem  to 
give  a  good  fit  to  the  data  points.  Let  R ,  represent  the  ratio  of  integrated  bed  stress,  then  for  V>0 

(4) 

where  a  was  found  to  be  403.  For  V>0,  appeared  to  behave  like 

Ri  = - l—z~  (5) 

b+c  ln(V) 

where  h  *  as  found  to  be  1.506  and  c  =  .1707.  While  this  last  relation  performs  well  for  the 
larger  values  of  V,  it  is  not  defined  for  no  ventilation  and  gives  /?,=1  for  a  finite  value  of  V.  For 
the  purposes  of  the  figure,  a  straight  line  has  been  drawn  between  the  point  at  V  =  .00025  and  V  = 
0. 


Even  in  a  turbulent  boundary  layer,  viscous  stresses  dominate  within  the  viscous  sub¬ 
layer.  This  means  that  the  bed  stress  can  be  approximated  as 


du 

dz 
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if  u  is  measured  sufficiently  close  to  the  bed.  It  is  generally  accepted  that  the  viscous  sublayer  is 
present  for  z  +  <  5  where  z*  is  the  dimensionless  vertical  coordinate  defined  by 


z 


+ 


u*z 


V 


(7) 


26 


and  u*  is  the  friction  velocity  vt/p .  The  lowest  elevation  for  which  velocity  was  measured  was 
at  0.01  cm  which  is  within  the  viscous  sublayer  for  all  phases.  It  is  therefore  possible  to  calculate 
a  second  estimate  of  bed  stress  using  (7).  This  process  has  been  carried  out  using  the  velocity 
record  at  0.01  cm  and  the  results  have  been  plotted  on  Figure  10  as  filled  boxes  and  show  good 
agreement  with  the  shear  stress  sensor  data 

4.  PERSISTENCE  OF  VENTILATION  EFFECT 

Re  Dependence 

In  order  to  determine  the  Re  dependence  of  the  effects  of  boundary  ventilation,  a  series 
of  tests  were  performed  in  which  the  flow  Re  was  varied.  Using  the  linear  relationship  between 
d0  and  um,  it  can  be  shown  that  (1)  is  equivalent  to 

Re  =  doa/4v.  (8) 

As  all  tests  were  performed  using  the  same  frequency,  reductions  in  Re  were  achieved  by  reduc- 
ing  d o  by  the  appropriate  factor.  The  range  over  which  Re  was  varied  was  determined  by  limits 
of  the  physical  apparatus.  The  largest  Re  (l.OxlO6)  represented  the  largest  stroke  which  could 
be  safely  obtained  from  the  hydraulic  ram  and  the  smallest  (l.OxlO5)  represented  the  lowest 
volume  of  water  which  could  be  accurately  pumped  through  the  permeable  bed  at  the  V  value  of 
lxlO-3. 

Figure  1 1  shows  the  results  of  these  tests.  At  each  Re  the  ensemble  averaged  mean 
stress  with  ventilation  has  been  integrated  over  one  oscillatory  half  cycle.  These  values  were  then 
normalized  by  the  integrated  bed  stress  for  the  unvcntilated  case  at  the  appropriate  Re.  These 
values  have  been  plotted  as  the  circles  and  squares  in  Figure  11.  The  bars  on  the  plots  represent 
95%  confidence  intervals  and  the  lines  are  linear  least  squares  fits  to  the  points.  As  can  be  seen  in 
the  figure  any  Re  dependence  is  weak  as  the  fit  to  the  data  are  essentially  flat.  Neither  slope  is 


significantly  different  from  0  at  the  95%  confidence  level  although  the  slope  for  suction  is 
significant  at  the  90%  confidence  level. 

Phase  Dependence 

Since  most  of  the  results  presented  in  this  paper  are  for  boundary  ventilation  which  is  in 
phase  with  the  free  stream  flow,  it  was  decided  to  test  the  sensitivity  of  the  flow  to  variations  in 
the  phase  relationship  between  the  two  flows.  Zero  phase  in  this  section  is  defined  as  when  the 
maximum  magnitude  of  the  boundary  ventilation  occurs  simultaneously  with  the  maximum  in  die 
free  stream.  Positive  phase  implies  that  that  maximum  ventilation  occurs  prior  to  the  maximum 
free  stream  velocity  and  negative  phase  therefore  implies  maximum  ventilation  follows  the  free 
stream  maximum.  The  range  of  phases  has  been  restricted  to  ±90°  with  separate  tests  for  suction 
and  injection. 

It  should  be  noted  that  any  phase  difference  between  the  two  flows  implies  that  each  half 
cycle  will  experience  periods  of  both  suction  and  injection.  For  that  reason,  it  was  determined 
that  presenting  results  in  terms  of  bed  stress  integrated  over  a  half  cycle  would  not  be  useful  and 
instead  results  are  presented  in  terms  of  maximum  bed  stress  normalized  by  maximum  unven¬ 
tilated  bed  stress.  These  values  are  presented  in  Figure  12  where  open  circles  represent  the 
results  for  suction  and  the  squares  are  the  values  for  injection.  The  error  bars  represent  95% 
confidence  limits  calculated  from  the  variance  of  the  maximum  stress  measurements. 

The  simplest  prediction  for  the  effect  of  phase  would  indicate  that  flow  history  has  no 
role  in  the  effect  of  boundary  ventilation.  Such  a  formulation  would  state  that  the  only  parameter 
of  importance  is  the  instantaneous  ventilation  parameter  V'.  This  parameter  can  be  easily  calcu¬ 
lated  as  the  cosine  of  the  phase  times  the  flow  ventilation  parameter  V.  Using  these  values  of  V\ 
the  expected  stress  reduction  or  enhancement  can  be  predicted  from  a  plot  similar  to  Figure  10 
which  is  based  on  stress  maxima  rather  than  integrated  stress.  This  procedure  was  used  in 
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constructing  the  lines  in  Figure  12. 

In  a  gross  sense  the  results  do  appear  to  follow  this  trend.  The  effect  of  ventilation  is  a 
maximum  around  a  phase  of  0°  and  drops  off  to  no  effect  (ratio  =  1)  at  ±90°.  The  effects  arc 
seen  not  to  be  a  critical  function  of  phase,  that  is  phase  differences  of  ±10°  do  not  radically 
change  the  effects  of  ventilation.  In  fact,  the  changes  in  bed  stress  over  this  range  are  not  statisti¬ 
cally  different  from  each  other.  However  for  larger  phase  differences,  the  drop  off  rate  is  clearly 
different  for  each  quadrant  and  therefore  warrants  inspection  in  all  quadrants.  The  most  straight 
forward  behavior  occurs  for  boundary  suction  with  positive  phase.  Here  the  effects  of  ventilation 
appear  to  follow  the  simple  cosine  response.  This  indicates  that  with  a  stabilized  boundary  layer 
as  exists  with  suction  the  stream  lines  equilibrate  rapidly  toward  their  undistu.bcd  state  as  the 
level  of  applied  suction  is  reduced.  While  not  behaving  as  smoothly  as  suction  with  positive 
phase,  injection  with  negative  phase  appears  to  also  follow  the  cosine  behavior.  It  is  important  to 
recall  the  sequence  in  the  development  of  the  boundary  layer  for  this  case.  In  particular  the  half 
cycle  starts  with  suction  occurring  in  the  boundary  which  switches  at  some  point  to  injection.  We 
have  already  seen  (above)  that  the  boundary  layer  will  quickly  respond  to  the  reduction  of  suction 
so  the  tendency  for  the  bed  stress  to  follow  the  cosine  response  implies  that,  when  applied  to  a 
stable  or  unventilated  boundary  layer,  the  response  time  for  shear  reduction  due  to  injection  is 
small  with  respect  to  the  time  scales  of  the  flow. 

The  behavior  of  injection  with  positive  phase  contrasts  strongly  with  the  above  descrip¬ 
tions.  In  fact  the  only  resemblance  that  the  effect  of  boundary  injection  with  positive  phase  has 
to  cosine  response  is  that  the  stress  ratio  docs  return  to  one  at  90°.  Figure  12  illustrates  that 
stress  reduction  due  to  injection  grows  at  phases  greater  than  0°  even  though  the  instantaneous 
ventilation  parameter  is  decreasing.  To  understand  this  behavior  it  is  important  once  again  to 
consider  the  time  history  of  the  boundary  layer.  In  this  situation,  the  boundary  layer  has  experi¬ 
enced  boundary  injection  for  the  entire  half  cycle  prior  to  the  occurrence  of  maximum  stress.  The 
maximum  injection  rate  occurs  prior  to  the  occurrence  of  maximum  stress  which  occurs  during  a 
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time  of  decreasing  injection.  With  this  scenario  in  mind,  it  appears  that  the  spreading  of  boun¬ 
dary  layer  streamlines  induced  by  boundary  injection  is  resistant  to  readjustment.  In  fact  the 
slight  tendency  for  increased  stress  reduction  with  phase  suggest  that  once  a  level  of  spreading  is 
developed,  it  can  be  maintained  by  smaller  ventilation  flows  than  originally  present.  In  contrast  to 
compressed  boundary  layer  which  appear  to  readjust  immediately  to  changing  conditions,  boun¬ 
dary  layer  spreading  appears  to  be  resistant  to  readjustment.  This  result  is  fully  consistent  with 
results  for  steady  flow  subject  to  changing  transpiration. 

This  interpretation  seems  to  be  validated  by  the  results  for  boundary  suction  with  nega¬ 
tive  phase.  In  this  quadrant,  the  half  cycle  commences  with  boundary  injection  leading  to  boun¬ 
dary  layer  spreading.  With  the  onset  of  suction,  the  stream  lines  are  indeed  drawn  closer  to  the 
bed  leading  to  greater  bed  stress  for  greater  suction,  but  the  level  of  bed  stress  enhancement  never 
manages  to  attain  the  level  it  does  working  on  a  boundary  layer  which  was  not  initially  subject  to 
injection. 

Sensitivity  to  Velocity  Asymmetry 

In  order  to  test  how  sensitive  the  effects  of  boundary  ventilation  are  to  the  shape  of  the 
free  stream  velocity,  a  series  of  experiments  were  performed  with  an  asymmetrical  velocity 
waveform,  a  plot  of  which  is  shown  in  Figure  1 3.  This  waveform  was  taken  from  a  time  series  of 
orbital  velocities  beneath  near  breaking  waves  in  the  ocean.  This  particular  wave  was  chosen  due 
to  it’s  strong  asymmetry,  seven  second  period  and  the  similarity  of  its  flow  parameters  to  the 
oscillatory  wave  forms  generally  employed  in  this  study.  The  orbital  displacement  distance  d o 
for  this  form  was  182  cm.  The  maximum  orbital  velocity  was  1 1 1  cm/s.  Using  (1)  this  give  a  Re 
of  l.OxlO6. 

A  series  of  runs  were  carried  out  to  test  the  affects  of  boundary  ventilation  on  this 
waveform.  In  these  tests  the  exact  same  wave  form  was  used  to  drive  the  ventilation  velocity  as 
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was  used  to  drive  free  stream  velocity  and  no  phase  difference  was  imposed.  It  should  be  under¬ 
stood  that  all  suction  results  come  from  the  peaked  high  velocity  portion  of  the  waveform  (crest) 
and  all  injection  results  come  from  the  flat  low  velocity  portion  of  the  waveform  (trough).  Once 
again,  the  integrated  bed  stress  for  each  sub  cycle,  normalized  by  the  unventilated  integrated  bed 
stress  for  the  appropriate  sub  cycle  is  plotted  as  the  open  symbols  in  Figure  1 3.  The  solid  curve  is 
the  same  one  plotted  in  Figure  10.  It  can  be  seen  from  this  figure  that  in  general  shape  has  an 
insignificant  affect  on  the  effects  of  boundary  suction.  Shape  does  seem  to  have  some  affect  on 
the  effect  of  injection,  particularly  for  larger  values  of  ’/.  This  can  be  understood  from  the  dis¬ 
cussion  in  the  previous  section.  It  can  be  seen  that  the  flat  profile  of  the  low  velocity  section 
leads  to  injection  approaching  its  maximum  value  early  in  the  sub  cycle  and  then  maintaining  that 
value  for  a  significant  portion  of  the  sub  cycle.  It  seems  likely  that  injection  in  this  case  would 
result  in  flow  streamlines  which  exhibit  a  higher  mean  separation  over  the  course  of  the  longer 
duration  sub  cycle  than  would  be  expected  with  a  more  sinusoidal  profile.  This  would  result  in 
greater  bed  stress  reduction  as  is  seen  in  the  figure. 

5.  DISCUSSION 

The  results  of  these  experiments  show  that  the  ventilated  oscillatory  boundary  layer 
exhibits  behavior  which  could  have  been  qualitatively  predicted  from  knowledge  of  the  manner 
in  which  transpiration  affects  steady  boundary  layers.  The  importance  of  these  effects  can  take  on 
a  larger  significance  in  the  oscillating  case  than  they  would  appear  to  in  the  steady  case.  This  can 
be  demonstrated  by  considering  the  change  in  the  boundary  layer  velocity  profile.  In  an  oscilla¬ 
tory  boundary  layer  where  the  velocities  in  one  half  cycle  are  just  the  negative  of  the  velocities  in 
the  second  half  cycle,  the  net  velocity  over  one  full  cycle  is  identically  zero.  However  in  a  ven¬ 
tilated  oscillatory  boundary  layer  where  suction  leads  to  higher  velocities  drawn  nearer  to  the  bed 
while  injection  causes  a  spreading  of  the  velocity  profile  with  lower  velocities  near  the  bed,  a  net 
velocity  results  in  the  boundary  layer.  The  sease  of  this  velocity  is  in  the  direction  of  the 
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oscillatory  How  which  is  subject  to  suction. 

A  mean  net  boundary  layer  velocity  has  been  determined  for  the  values  of  V7  for  which 
velocity  profiles  where  measured.  This  was  done  by  integrating  the  mean  velocity  profiles  over  a 
height  of  2  So  (3.6  cm)  where  So  was  defined  earlier.  This  vertically  averaged  mean  velocity  was 
then  integrated  over  an  entire  cycle  to  arrive  at  the  mean  vertically  averaged  net  velocity  herein 
called  the  ventilation  current  <v>v.  These  values  have  been  plotted  in  Figure  14.  The  line  in 
Figure  14  is  a  least  squares  fit  for  the  relation  <u>v/um=qVn.  The  best  fit  was  found  to  be 

<u>vlum-  4.64x |  V7 1 °'6 .  (9) 

In  this  relation,  the  sign  of  <u  >v  is  opposite  the  sign  of  V. 

While  the  bed  experiences  no  net  stress  in  a  purely  oscillatory  boundary  layer,  there  is  a 
net  bed  stress  with  ventilation.  This  net  stress  is  a  consequence  of  the  bed  stress  reduction  due  to 
injection  and  the  bed  stress  enhancement  due  to  suction.  When  the  two  half  cycles  are  averaged 
together  they  result  in  a  net  bed  stress.  Figure  IS  shows  the  net  bed  stress  due  to  a  ventilated 
oscillatory  boundary  layer  as  a  function  of  V.  Here  the  results  of  Figure  9  have  been  used  to  con¬ 
struct  this  plot.  The  points  in  Figure  14  represent  the  integrated  difference  in  the  bed  stress 
between  the  bed  stress  with  suction  and  the  bed  stress  with  injection  normalized  by  twice  the 
integrated  unventilated  bed  stress.  These  points  therefor  represent  the  mean  net  stress  over  an 
entire  period  as  a  percentage  of  the  mean  gross  stress  for  the  unventilated  case.  The  solid  line  in 
the  figure  is  straight  line  interpolation  between  the  points.  The  dashed  line  is  a  linear  fit  to  the 
first  three  values  for  ventilation  which  was  forced  to  go  through  the  origin.  It  appears  that  for 
small  ventilation  (V<0.001),  the  net  stress  can  be  approximated  by  the  linear  relation 

<t>v/<|tC)|>  =  -420xV7.  (10) 

Notice  that  this  relation  suggests  that  a  1%  net  stress  is  obtained  with  a  ventilation  parameter  as 
small  as  2.4xl0~5!  It  should  again  be  observed  that  the  direction  of  this  net  stress  is  in  the 
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direction  of  the  oscillating  flow  experiencing  suction. 

It  is  worthwhile  to  consider  the  importance  these  results  may  hold  for  geophysical  situa¬ 
tions.  Reid  &  Kajiura  [1957]  as  well  as  others  have  shown  that  a  gravity  wave  traveling  through 
a  fluid  medium  over  a  permeable  matrix  will  induce  boundary  ventilation.  It  can  be  shown  that 
for  geophysically  reasonable  permeabilities,  the  ventilation  parameter  V  will  not  exceed  1x10 
and  that  phase  differences  between  u  and  w  will  be  O  (5°)  or  less.  This  means  that  in  general 
wave  driven  ventilated  boundary  layers  will  produce  a  mean  boundary  current  in  the  direction  of 
wave  advance  and  with  a  magnitude  which  can  be  determined  through  a  knowledge  of  V.  A  net 
stress  in  the  direction  of  wave  advance  will  be  induced,  with  a  magnitude  which  can  be  predicted 
by  (10).  These  phenomena  will  be  Re  independent  over  an  order  of  magnitude  of  Re  as  shown  in 
Figure  11.  Over  the  expected  range  of  phase  separation  between  the  wave  potential  flow  and 
boundary  ventilation,  phase  would  not  be  expected  to  significantly  alter  these  results.  While 
gravity  waves  are  not  expected  to  exhibit  a  sinusoidal  profile  and  will  develop  inherent  asym¬ 
metries,  the  results  of  Figure  13  indicate  that  for  low  values  of  V  the  flow  will  still  exhibit  the 
above  discussed  behavior. 

This  study  shows  that  ventilated  oscillatory  boundary  layers  will  develop  net  boundary 
velocity  as  well  as  net  bed  stress  even  in  the  absence  of  any  other  asymmetry.  This  result  has 
obvious  implications  for  studies  of  transport  in  such  flows.  The  implications  of  the  asymmetry  in 
turbulence  is  not  so  obvious  to  predict.  Figure  7  shows  quite  clearly  that  the  majority  of  the  tur¬ 
bulence  throughout  the  cycle,  originates  from  the  injection  half  cycle.  Yet  how  this  would  affect 
the  transport  of  some  quantity  is  far  from  clear.  If  transport  is  approximated  by  the  product  of 
suspension  and  local  velocity  and  the  level  of  suspension  is  proportional  to  the  instantaneous  tur¬ 
bulence  levels,  the  argument  can  be  made  that  transport  would  once  again  be  in  the  direction  of 
flow  during  suction.  This  can  be  understood  by  examining  Figure  16.  This  figure  gives  a  time 
history  of  the  vertically  averaged  turbulence  levels  present  in  Figure  7.  This  figure  which  imparts 
no  information  about  history  clearly  shows  that  the  maximum  vertically  averaged  turbulence 


levels  occur  during  the  suction  half  cycle.  By  the  above  assumptions,  transport  due  to  turbulent 
suspension  would  be  in  the  direction  of  flow  during  suction  even  in  the  absence  of  other  asym¬ 
metries.  As  suspended  transport  under  waves  is  a  more  complicated  phenomenon  than  simple 
instantaneous  averages  the  above  scenario  may  be  wholly  inaccurate.  However  it  does  help  to 
demonstrate  the  possible  ramifications  these  findings  have  for  suspended  transport. 

6.  CONCLUSIONS 

A  laboratory  experiment  was  performed  to  test  the  affect  of  boundary  ventilation  on  the 
mean  velocities,  turbulent  flow  characteristics  and  bed  stresses  in  an  oscillatory  boundary  layer. 
The  results  of  these  experiments  lead  to  the  following  conclusions: 

1.  Boundary  suction  causes  flow  streamlines  to  be  pulled  toward  the  bed,  leading  to  higher 
velocities  near  the  bed,  while  injection  lead  to  streamline  spreading  and  reduced  near  bed  veloci¬ 
ties.  This  effect  results  in  a  mean  near  bed  flow  in  the  direction  of  the  flow  experiencing  suction. 

2.  The  enhanced  velocity  shear  near  the  bed  during  suction  resulted  in  increased  bed  stress 
while  the  reduced  shear  during  injection  results  in  diminished  bed  stress. 

3.  The  above  affects  are  shown  to  be  strong  functions  of  the  ventilation  parameter  V.  The 
bed  stress  is  shown  to  be  independent  of  Re  for  the  range  of  Re  tested.  It  is  also  shown  that 
phase  dependence  is  small  for  small  phase  differences  between  u  and  w.  The  effects  of  suction 
were  insensitive  to  waveform  shape  while  the  effects  of  injection  increased  with  increasing  length 
of  the  injection  sub-cycle.  This  tendency  was  negligible  for  small  values  of  V. 

4.  Boundary  ventilation  leads  to  a  strong  asymmetry  in  flow  turbulence.  Enhanced  tur¬ 
bulent  fluctuations  are  maintained  in  a  thin  layer  near  the  bed  with  suction  and  turbulent  diffusion 


is  reduced.  Away  from  the  bed  turbulence  levels  are  enhanced  by  injection  leading  to  greater  tur¬ 
bulent  d.ffusion  initiated  during  injection.  The  result  of  these  effects  is  enhanced  vertically  aver¬ 
aged  turbulence  levels  during  suc'ion. 
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Ensemble  averaged  time  series  of  (a)  centerline  velocity,  (b)  centerline 
acceleration,  and  (c)  bed  stress  for  the  unventilated  case.  Dashed  lines  in 
(c)  represent  ±  one  standard  deviation  in  absolute  terms.  Solid  circles  are 
bed  stress  estimates  from  log  layer  fit  using  ks  calculated  from  the  first  half 
cycle.  Open  circles  represent  bed  stress  estimates  using  ks  calculated  from 
2nd  half  cycle. 


CONTOURS  OFV  w'2  FOR  w/u  =  -0.001,  (l80°-0°  =  INJECTION) 


right  correspond  to  measurement  locations. 


ure  9:  Ensemble  averaged  bed  stress  for  various  values  of  V.  Solid  line  is  bed  stress  for  unventilated  case.  Brackets 

represent  957c  confidence  level  of  peak  stress. 


VENTILATION  PARAMETER  (w/u) 

Time  averaged  ventilated  bed  stress  vs  injection  parameter  V.  Values  are 
normalized  by  time  averaged  unventilated  bed  stress.  Averages  are  for  first 
half  cycle.  Brackets  represent  95%  confidence  interval  based  on  repeated 
estimates  of  ensemble  averaged  quantity.  Variance  is  assumed  to  be  a  con¬ 
stant  fraction  of  measurement.  Circles  are  for  VcO  and  boxes  are  for  V/>0. 
Solid  boxes  represent  estimates  based  on  du/dz  over  the  bottom  0.01  cm. 
Solid  line  is  from  equation  (4)  and  equation  (5). 


Phase  of  w  relative  to  u  (<t>) 

Maximum  ventilated  bed  stress  values  normalized  by  maximum  unven¬ 
tilated  bed  stress  as  a  function  of  ventilation  phase  <J>.  Open  circles  are  for 
suction,  open  boxes  are  for  injection.  Solid  lines  are  expected  response 
assuming  instantaneous  response.  Brackets  represent  95%  confidence 
interval  based  on  variance  of  ensemble  averaged  estimate. 
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VENTILATION  PARAMETER  (w/u) 

Figure  13:  Time  averaged  ventilated  bed  stress  vs  injection  parameter  V  for  asym¬ 
metric  wave  form  (symbols).  The  line  represents  results  for  a  symmetrical 
wave  form  (from  Figure  10).  Values  are  normalized  by  time  averaged 
unventilated  bed  stress.  Averages  are  over  one  half  cycle.  The  peaked  half 
cycle  for  V<()  and  flat  half  cycle  for  V>0.  Brackets  represent  95 9r 
confidence  interval  based  on  repeated  estimates  of  ensemble  averaged 
quantity.  Variance  is  assumed  to  be  a  constant  fraction  of  measurement. 
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Figure  14:  Net  boundary  layer  velocity  averaged  over  a  distance  250  (ventilation  current)  vs  absolute  value  of  V.  Net  velocity 
is  presented  as  a  percentage  of  Line  is  equation  (9).  Sign  of  net  velocity  is  opposite  the  sign  of  V. 
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Figure  15:  Average  net  ventilated  bed  stress  as  a  percentage  of  average  gross  unven 
tdated  bed  stress  vs  V  Dashed  line  is  linear  fit  to  points  for  VSlxIfP1 
Brackets  are  95 %  confidence  interval. 


Figure  16:  Vertically  averaged  turbulence  C*wu  )  level  as  a  function  of  wave  phase  0.  Average  is  over  8  cm  and  V'=lxlO 
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CHAPTER  III 


WAVE  DRIVEN  VENTILATION  IN  SAND  BEDS 

1.  INTRODUCTION 

Any  beach  which  is  composed  of  sediments  of  a  non  local  origin  stands  as  an  incontro¬ 
vertible  proof  of  the  existence  of  a  shoreward  directed  force  opposing  the  constant  down  slope 
and  hence  seaward  pull  of  gravity.  The  source  of  this  forcing  has  long  been  recognized  to  be  the 
waves  which  shoal  and  dissipate  on  these  beaches  but  the  mechanism  through  which  these  waves 
generate  this  forcing  is  the  subject  of  much  discussion.  The  role  that  wave  driven  ventilation  may 
play  in  this  process  is  investigated  in  this  study. 

A  surface  gravity  wave  traveling  in  a  fluid  overlying  a  permeable  bed  has  been  predicted 
to  induce  an  oscillating  flow  in  that  bed  [Putnam,  1949;  Reid  &  Kajiura,  1957],  This 
phenomenon  was  originally  of  interest  to  coastal  engineers  and  scientists  because  it  was  proposed 
as  a  mechanism  for  the  dissipation  of  wave  energy  prior  to  breaking.  More  recently  Conley  and 
Inman  [1992]  labeled  this  type  of  phase  locked  oscillatory  flow  in  the  bed  "ventilation",  and  sug¬ 
gested  that  boundary  layer  phenomenon  associated  with  ventilation  could  explain  an  asymmetry' 
which  they  observed  in  the  development  of  fluid  granular  boundary  layer  under  waves.  This 
work  reports  on  field  observations  of  wave  driven  ventilation,  compares  the  observations  to 
theory  and  discusses  the  importance  of  ventilation  in  nearshore  processes. 
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2.  VENTILATION  THEORY 


In  the  most  basic  treatment  of  this  problem  we  consider  a  fluid  of  depth  h  overlying  a 
bed  of  thickness  b  (Figure  1)  [e.g.  Putnam,  1949).  The  bed  is  assumed  to  be  homogeneous  with 
permeability  K  and  to  overlay  an  impermeable  substrata.  Flow  within  the  bed  is  assumed  to  fol¬ 
low  Darcy’s  law  so  that 

-K  fyb  -K  dPb 

u  - - —  w  = - (1) 

g  dx  | X  dz 


where  ub  and  Wf,  are  the  vertical  and  horizontal  flow  in  the  bed  respectively,  p b  is  pressure  in  the 
bed  and  |i  is  the  viscosity  of  the  pore  fluid.  Now  substituting  equation  (1)  into  the  continuity  rela¬ 
tion  gives 


d2pb  &Pb 

dx 2  az2 

Equation  2  must  be  solved  subject  to  the  following  boundary  conditions 


(2) 


pb=p  @  z  =0 


(3) 


and 


wb  =  0  @  z  =  -b. 


(4) 


Here  p  is  the  pressure  in  the  fluid.  From  linear  wave  theory  it  is  known  that 


p=pag 


cosh ( kz ) 
cosh (kh ) 


cos  ( at-kx ) 


(5) 


where  a  is  the  wave  amplitude,  k  =  2: l/L  is  the  wave  number,  a  =  2n/T  is  the  wave  radian  fre¬ 
quency  and  L  and  T  are  the  wavelength  and  period  respectively.  The  solution  for  (2)  which 
satisfies  (3)  and  (4)  is 


Pb  = 


Pan 

cosh ( kh ) 


cosh  ffc  (z  +b )] 
cosh (kb ) 


cos  (at-kx). 


(6) 


From  (1 )  and  (6)  we  can  see  that  flow  is  induced  in  the  bed  by  the  pressure  field  associated  with 


the  wave  and  that  the  vertical  ventilation  will  be  of  the  form 


-K  agk  sinh[k(z+b)]  . 

wb  = - T7IIT - cos(ot-kx). 

v  cosh(kh)  cosh  (kb) 


(7) 


This  flow  is  an  oscillating  flow  into  and  out  of  the  bed  and  is  phase  locked  to  the  wave  pressure 
such  that  maximum  flow  into  the  bed  (wb<0)  occurs  under  the  crest  of  the  wave  and  maximum 
flow  out  of  the  bed  occurs  under  the  trough  of  the  wave  (Figure  2).  A  similar  relationship  exists 
between  vertical  ventilation  wb,  and  the  horizontal  wave  orbital  motion  in  the  fluid  u,  where 

k  cosh(kz) 


u~ag - --■■■■--cos  (Gt-kx). 

ct  cosh(kh) 


(8) 


Comparison  of  (7)  and  (8)  shows  that  for  linear  progressive  waves  maximum  ventilation  into  the 
bed  occurs  simultaneously  with  maximum  horizontal  orbital  motion  under  the  crest  (onshore 
orbital  motion)  and  maximum  ventilation  out  of  the  bed  occurs  simultaneously  with  maximum 
offshore  directed  orbital  motion.  Because  orbital  velocities  and  ventilation  are  in  phase,  the  ratio 
Wb  /  u  is  a  constant.  At  the  fluid  bed  boundary  (z  =0),  this  ratio  is 


—  =  -—or  tanh  (kb).  (9) 

w  v 

Here  v  =  p./p  is  the  kinematic  viscosity.  Inspection  of  (9)  suggests  the  relative  magnitudes  of 
Wb  and  u.  In  particular,  for  common  wave  frequencies  in  the  nearshore,  a  tanh  (kb)  will  be 
assumed  to  be  about  1 .  Darcys  law  is  generally  assumed  valid  where  the  grain  Reynolds  number, 
Re^  =  wbD  N.  is  O  (1)  or  less  where  D  is  the  grain  diameter.  It  is  reasonable  to  make  «=102 
cm/s,  then  with  v=.01cm2/s,  (9)  gives  wb  =  10 4xK.  Using  this  in  the  Reg  limit  gives 
KD  =  10-6cm3.  Krumbcin  and  Monk  [1942]  have  shown  that  the  permeability  in  sands  is 
proportional  to  the  square  of  the  mean  grain  diameter.  Explicitly, 

K  =  7.6xl0_4£>2e^-1'3'a*)  (10) 

where  D  is  in  cm  and  O0  is  the  deviation  measure  in  phi  units.  The  result  is  permeability  in  cm2. 
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Approximating  this  relation  as  K  =  lO^D  2  the  Re^  limit  is  now  ax  -  10_2cm3 .  In  order  to 
determine  the  largest  value  (9)  could  take,  we  shall  let  D^ax  =  10-1cm2  which  when  combined 
with  (10)  and  (9)  tells  us  that  wb/u  is  <9(1(T3)  or  less.  Notice  that  the  above  discussion  indi¬ 
cates  that  (1)  will  be  valid  for  grain  sizes  on  the  order  of  mm  and  smaller.  This  is  consistent  with 
the  findings  of  Gu  and  Wang  [1991]  and  others. 

3.  EXPERIMENTAL  PROCEDURE  AND  RESULTS 

Using  the  order  of  magnitude  arguments  mentioned  above  for  K  and  v,  it  can  be  seen 
from  (1)  that  the  vertical  pressure  gradient  in  the  bed  is  expected  to  be  O  (103)  times  the  ventila¬ 
tion.  It  was  therefore  decided  that  ventilation  in  field  conditions  could  best  be  determined  from 
measurements  of  the  vertical  pressure  gradient  in  the  bed.  The  procedure  was  to  bury  a  differen¬ 
tial  pressure  sensor  in  a  sand  bed  with  the  pressure  ports  separated  in  the  vertical  direction  (Fig¬ 
ure  3).  Such  an  arrangement  results  in  a  steady  pressure  differential  equal  to  the  hydrostatic  head 
between  the  two  ports  as  well  as  a  fluctuating  component  which  is  directly  proportional  to  the 
ventilation.  Pressure  gradients  due  to  the  slope  of  the  wave  are  effectively  filtered  out  as  there  is 
no  separation  of  the  ports  in  the  crosshore  direction. 

The  experiment  was  performed  on  20  and  21  August  1992  near  the  Scripps  Institution  of 
Oceanography  pier.  A  Setra  model  228-1  differential  pressure  sensor  with  a  port  separation  of  10 
cm  was  buried  in  the  sand  bed.  The  top  port  was  placed  6  cm  below  the  sand  water  interface.  A 
Marsh-Mcbimey  model  512  two  component  electromagnetic  (em)  current  meter  and  a  Setra 
model  280  absolute  pressure  sensor  were  placed  in  the  water  column  directly  above  the  buried 
sensor.  Data  was  collected  at  8  Hz  over  a  two  and  one  half  hour  run.  Five  runs  were  made  span¬ 
ning  an  entire  tidal  cycle  with  mean  water  depth  ranging  from  290  to  340  cm.  The  local 
significant  wave  height  ranged  from  50  -  66  cm  at  a  peak  wave  period  of  6.5  s.  Sand  in  the  top  of 
the  bed  had  a  mean  grain  diameter  of  180  microns  with  a  0$  of  0.55.  This  gives  a  permeability 
of  12  Darcys  (10xl(T8cm2)  using  (10).  Cores  were  taken  at  the  instrument  burial  site  during  the 
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experiment  and  were  later  subject  to  a  hydrostatic  head  in  order  to  determine  the  permeability  by 
direct  measurement.  This  technique  gave  a  permeability  of  1 1  Darcys  which  is  the  value  used  in 
this  work. 

Figure  4  shows  a  typical  section  of  time  series.  The  top  trace  represents  sea  surface 
elevation  as  measured  by  the  absolute  pressure  sensor,  the  middle  trace  is  crosshore  orbital 
motion  ( u )  46  cm  above  the  bed  as  measured  by  the  em  current  meter  and  the  bottom  trace  is  ven¬ 
tilation  as  calculated  from  the  differential  pressure  sensor  measurements.  It  is  seen  that  ventila¬ 
tion  is  present  and  behaves  qualitatively  as  predicted  by  (7).  In  particular  maximum  ventilation 
into  the  bed  occurs  under  the  crest  of  the  wave  concurrently  with  maximum  onshore  orbital 
motion.  Likewise  maximum  flow  out  of  the  bed  occurs  under  the  trough  of  the  wave  concurrently 
with  maximum  offshore  orbital  motion.  It  can  also  be  seen  that  the  ratio  of  |  wb/u  |  is  O(l0~5). 

If  equation  (7)  is  correct  then  it  should  be  possible  to  use  (5)  to  predict  the  ventilation  at 
a  depth  -zb  from  knowledge  of  the  pressure  at  a  height  above  the  bed  equal  to  zp.  The  relation 
would  be 


Wfa(Zft)  =  ■ 


K,  sinh[k(b-zb )] 

- fa - 

p.  cosh  ( kzp)cosh  {kb) 


P{zp). 


(11) 


As  the  processes  discussed  arc  linear,  it  should  be  possible  to  construct  the  power  spectrum  for 
wb  from  p  using  (11).  The  only  parameters  in  (11)  which  are  not  fixed  by  either  the  wave  fre¬ 
quency  or  fluid  properties  are  the  bed  permeability  K  and  the  bed  thickness  b.  However  as  K  is 
used  in  our  determination  of  wb,  the  only  previously  unspecified  parameter  in  this  problem  is  b. 

Figure  *  g'vp.s  »h<>  results  of  such  a  construction.  The  solid  line  in  this  figure  is  the  spec¬ 
tra  of  wb  as  calculated  from  the  differential  pressure  sensor  response.  The  dotted  line  is  the 
predicted  spectra  for  wb  as  calculated  from  the  absolute  pressure  spectra  and  a  bed  thickness  of 
250  cm.  This  thickness  represents  the  minimum  average  sand  bed  thickness  known  to  occur  in 
this  area.  It  can  be  seen  from  the  figure  that  using  this  depth  results  in  a  significant  under 
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prediction  of  ventilation.  It  was  found  that  to  best  match  the  peak  of  the  ventilation  spectra,  a  bed 
thickness  of  580  cm  was  required.  The  predicted  spectra  using  a  bed  thickness  of  580  cm  is 
shown  in  the  figure  as  the  dashed  line  and  it  can  be  seen  that  (7)  predicts  the  behavior  of  ventila¬ 
tion  at  least  over  the  energetic  part  of  the  spectrum.  This  thickness  (580  cm)  worked  equally  as 
well  for  all  five  runs. 

Knowledge  of  the  local  sediment  cover  indicates  that  580  cm  is  an  unrealistic  estimate 
of  the  sand  bed  thickness  but  consideration  of  the  local  beach  stratigraphy  may  give  an  explana¬ 
tion  of  where  this  value  comes  from.  It  is  know  that  the  sand  in  this  bed  overlies  a  cobble  layer 
about  50  cm  thick  (Figure  6).  If  we  assume  that  the  sand  bed  in  this  region  was  at  its  thickest 
summer  time  configuration,  a  sand  thickness  of  350  cm  is  not  unreasonable.  Add  50  cm  of  cob¬ 
ble  and  we  have  a  permeable  bed  of  400  cm.  This  is  still  short  of  the  required  580  cm  but  it  is 
suggested  that  this  can  be  accounted  for  by  considering  that  the  permeability  in  the  cobble  layer 
is  much  higher  than  that  in  the  sand  layer.  The  derivation  of  (7)  assumed  a  homogeneous  perme¬ 
able  bed  of  constant  permeability  K.  It  is  proposed  that  580  cm  represents  an  effective  bed  thick¬ 
ness  be.  That  is  that  the  two  layer  bed  of  sand  and  cobble  "effectively"  acts  likes  a  homogeneous 
sand  bed  of  permeability  K  with  a  thickness  of  about  580  cm. 

4.  DISCUSSION 

It  has  been  shown  that  ventilation  is  present  in  permeable  beds  beneath  waves  and  that 
it’s  behavior  can  be  well  predicted  by  simple  linear  theory  given  that  the  local  effective  bed  depth 
is  known.  It  would  now  be  useful  to  examine  what  significance  this  holds  for  sediment  transport 
and  other  nearshore  processes. 

Conley  and  Inman  [submitted!  have  shown  that,  similarly  to  steady  flow,  small  boun¬ 
dary  normal  velocities  can  strongly  affect  the  bed  stress  and  turbulent  boundary  layer  characteris¬ 
tics  of  an  oscillating  flow.  They  defined  a  ventilated  oscillatory  boundary  layer  as  the  boundary 
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layer  arising  from  a  purely  sinusoidal  bed  parallel  flow  subject  to  a  much  smaller  bed  normal  flow 
which  was  also  varying  sinusoidally  with  the  same  period  and  phase  as  the  larger  flow.  This  is 
just  the  situation  suggested  by  (7).  It  was  shown  that  the  bed  stress  arising  from  such  a  flow  when 
integrated  over  a  full  cycle  would  result  in  a  net  bed  stress  in  the  direction  of  the  flow  during  the 
half  cycle  experiencing  flow  into  the  bed.  This  concept  can  be  better  understood  by  examining 
the  schematic  of  wave  driven  ventilation  in  Figure  7.  Row  into  the  bed  under  the  crest  of  the 
wave  leads  to  an  enhancement  of  the  bed  stress  while  flow  out  of  the  bed  leads  to  a  reduction  in 
the  bed  stress.  The  effect  of  this  enhancement  and  reduction  averaged  over  one  wave  cycle  is  to 
produce  a  net  bed  stress  in  the  direction  of  wave  advance  or  an  onshore  directed  net  stress.  Con¬ 
ley  and  Inman  [submitted]  also  demonstrated  that  the  magnitude  of  the  normalized  net  stress  was 
a  function  of  w/u  and  that  for  small  values  of  w/u  (w/u<  10~3)  this  relationship  could  be 
approximated  by  the  empirical  linear  relationship 


<T>V 

<hol> 


= -420x 


(12) 


where  T  is  bed  stress,  Tq  is  the  bed  stress  without  ventilation,  <  >  implies  the  time  average  over 
one  cycle  and  the  subscript  v  indicates  a  ventilation  related  process.  Notice  that  <x>v  is  the  net 
bed  stress.  Using  (9)  and  (12)  we  have 

<T>V  —  420~cj  tunh  (kb)<  ] Tq  |  >  (13) 


A  kinematic  discussion  of  sediment  transport  requires  knowledge  of  the  time  history  of 
stress  and  velocity.  This  is  clear  when  one  considers  that  even  the  most  basic  sediment  transpon 
models  assume  that  transport  is  proportional  to  the  product  of  velocity  and  stress.  However  from 
a  force  balance  point  of  view,  net  stresses  are  very  important  in  nearshore  processes.  As  dis¬ 
cussed  earlier,  the  argument  here  states  that  in  a  situation  where  a  sloping  bed  of  sediment  is  con¬ 
stantly  agitated,  a  steady  force  must  be  exerted  in  the  up-slope  direction  to  counter  balance  the 
down-slope  force  due  to  gravity.  The  net  stress  due  to  boundary  ventilation  is  one  possible 
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source  for  just  such  a  force. 

Several  other  mechanisms  have  been  proposed  as  sources  for  such  a  force  and  it  is  of 
interest  here  to  make  a  comparison  of  the  importance  of  some  of  these  mechanisms  relative  to  the 
boundary  ventilation  mechanism.  The  focus  of  this  discussion  will  be  on  the  two  most  com¬ 
monly  mentioned  mechanisms;  namely  the  net  stress  arising  from  velocity  asymmetries  in  the 
nonlinear  shoaling  wave  form  [e.g.  Grant,  1943;  Inman  and  Bagnold,  1963]  and  the  net  stress 
arising  from  a  sinusoidal  oscillatory  flow  which  has  a  steady  boundary  current  imposed  on  it  as 
would  be  predicted  by  bottom  wind  theory  [Longuett-Higgins,  1953]. 

In  order  to  arrive  at  a  net  stress,  a  quadratic  stress  relation  is  assumed  where  the 
coefficient  of  friction  Cf  is  assumed  constant  through  out  the  entire  cycle.  This  is  represented  as 

x(0  =  p~«(0l«(0|.  04) 


where  the  period  —tp<t<tp  represents  the  portion  of  the  wave  period  T  over  which  velocity  u  is 
positive. 

To  represent  a  shoaling  gravity  wave,  we  use  the  third  order  Stokes  wave  as  expressed  in 
Koh  and  Le  Mehaute  [  1966].  The  wave  orbital  velocity  at  the  bed  is 

o  3 

«(/)  =  —  Y.  nsncos[n(kx-Gt)}  (16) 

kn  =  1 

where  sn  are  constants  determined  from  a,  k  and  h  and  listed  in  Koh  and  Lc  Mehaute  [1966].  It 
should  be  remembered  that  the  ratio  of  snlsn_ j  is  O  (e)  where  Z  —  ak  <t:l.  The  technique  used 
to  determine  tp  was  to  set  (16)  equal  to  0  and  solve  to  order  E2.  This  gives 


« 
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f  lsinhbkh  1 

sinh3  kh  ] 

{  V  9e2  2 

*  I 

Substituting  (16)  into  (15)  results  in 

<t>,  =  f  ^^|5o(4<,-rH4i  ^Jintoo/,)  j 

^0  =  5?  +4^2  +9^3 
£)  =45152  +  1252^3 

=51+65)53 

where :  £,3  =45)52 

£4  =452  +65)53 
^5  =  125253 

^6  =  M 


(17) 


(18) 


Here  tp  is  calculated  as  described  above. 

For  the  net  bed  stress  due  to  bottom  wind,  we  model  the  near  bed  velocity  as 


u(t)  =  umcos  (kx-Gt)  +  <u  >b 


(19) 


The  first  term  on  the  right  is  the  near  bed  velocity  from  linear  wave  theory  (8)  where 
um  =  ao/ sinh(kh)  is  the  maximum  near  bed  orbital  velocity.  The  second  term,  <u  >b ,  is  the 
steady  streaming  velocity  just  outside  the  boundary  layer,  the  "bottom  wind"  and  is  known  to 
consist  of  a  Eularian  and  a  Lagrangian  component  [Batchelor,  1985].  As  we  arc  discussing  a 
force  balance  at  a  fixed  point,  the  component  of  interest  is  the  Eularian  term  which  is 
(3/4)a2kcs!s\nh2(kh).  This  result  comes  from  calculations  for  a  laminar  boundary  layer  and 
laboratory  experiments  indicate  that  the  value  for  turbulent  boundary  layers  is  less.  Bijkcr  ct  al 
[1974]  report  that  <u>b  for  turbulent  flows  is  about  0.4  times  the  prediction  from  laminar 
theory.  Therefore,  for  this  work  the  value  used  will  be 


<u  >h  =  0.3 


a2ka 
sinh2  (kh) 


m 


m 


m 


m 


m 


(20) 
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The  equation  for  net  stress  due  to  bottom  wind  can  be  derived  using  (20),  (19)  and  (15)  as  well  as 
the  fact  that  tp  =  (l/a)cos-1(_<“  >;,/«„),  giving 


<T>b  = 


u2 

p 

nv 


1.+0. 


sin  ’(0.3C>+0.9!;Vl.-0.9C2 


(21) 


where : 


r  - 

19  sinh(kh). 


The  final  calculation  involves  deriving  an  expression  for  the  net  stress  due  to  ventilation. 
Once  again  we  assume  a  quadratic  stress  relationship  as  in  (14)  and  near  bed  velocity  as  in  (8). 
Fvaluating  the  mean  gross  stress  (<  |  to  |  >)  term  in  (13),  we  obtain  the  expression 

<T>V  =420— tanh(kb)^-CfU„o  (22) 

v  4  J 

It  is  helpful  to  normalize  (22)  in  order  to  separate  the  effects  directly  related  to  wave  intensity  and 
those  related  to  the  ventilation  process.  This  is  >  omplished  by  normalizing  by  the  peak  stress 
expected  without  ventilation  Tmax.  Inspection  of  (8)  and  (14)  show  that  xmax  =  Cfpu„/ 2.  The 
normalized  version  of  (22)  becomes 

<t>v  K 

- =  420—c  tanh  (kb)  (23) 

^max  2v 

As  can  be  seen  from  (23)  the  net  stress  due  »o  ventilation  is  approximately  a  linear  function  of 
frequency  and  permeability. 

Figure  8  shows  a  comparison  of  the  net  stress  due  to  these  three  generating  mechanisms 
as  a  function  of  grain  size  for  a  60  cm  high  wave  traveling  in  300  cm  deep  water.  The  net  stress 
due  to  the  ventilation  is  calculated  from  (23)  and  (10)  directly  while  the  net  stress  due  to  the 
Stokes  wave  asymmetry  and  bottom  wind  are  calculated  from  similarly  normalized  forms  of  (20) 
and  (21 )  respectively.  Figure  8a  shows  the  relative  strengths  of  net  stress  mechanisms  f.,r  wave 
parameters  similar  to  those  present  during  the  data  run  used  to  construct  Figure  5.  It  is  clear  from 
this  figure  that  the  effects  of  ventilation  are  strongly  grain  size  dependent.  For  fine  sand  the  net 
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bed  tress  due  to  the  ventilation  mechanism  appears  to  be  negligible  relative  to  the  oilier  two 
mechanisms.  For  medium  grain  sands  around  500  microns,  the  net  stress  due  to  each  of  the  three 
mechanisms  is  significant  and  for  coarse  grain  sands  and  pebbles  1  mm  in  diameter  or  greater,  net 
stress  due  to  ventilation  clearly  dominates.  The  affect  of  increasing  period  is  shown  in  Figures  8b 
and  8c.  It  is  ^een  that  ventilation  decreases  for  increasing  period  but  is  still  a  significant  factor  for 
coarse  grain  sediments. 

It  is  illuminating  to  consider  what  is  happening  when  the  grain  size  or  wave  frequency 
are  altered.  Comparing  (23)  and  (9),  we  see  that  changes  in  <x>v  come  from  changes  in  w/u. 
While  linear  theory  indicates  that  changing  the  grain  size  (permeability)  and  frequency  are  the 
only  means  to  change  w/u,  any  other  process  which  increases  this  ratio  will  result  in  increased 
importance  of  this  phenomenon.  One  such  mechanism  could  be  the  stressing  of  the  sediment  bed 
by  the  wave  orbital  motion.  The  bed  dilates  as  it  is  stressed  by  the  orbital  motion,  and  this  dila- 
‘.;on  leads  to  increased  pore  size.  This  increased  pore  size  results  in  increased  permeability  and 
therefore  increased  ventilation.  Maximum  stressing  of  the  bed  occurs  simultaneously  with  max¬ 
imum  ventilation.  Any  experiment  such  as  this  one  where  permeability  is  measured  indepen¬ 
dently,  would  not  detect  this  effect. 

This  limited  study  indicates  that  for  fine  sand  the  net  stress  due  ic  ven  ilation  is  less 
significant  than  that  due  to  bottom  wind  and  wave  form.  However,  this  does  not  mean  that  the 
affects  of  ventilation  are  negligible  in  crosshore  transport.  Conley  and  Inman  | submitted]  nave 
shown  that  another  of  the  effects  of  ventilation  is  a  disruption  in  the  development  of  turbulence  in 
the  oscillatory  boundarv  layer.  Their  results  indicate  that  during  suction  (flow  into  the  bed),  tur¬ 
bulence  is  drawn  closer  tr  the  bed  and  maintained  in  a  compact  layer  near  the  bed.  During  injec¬ 
tion  (flow  out  of  the  bed)  turbulence  extends  much  higher  into  the  fluid  bui  is  pushed  away  from 
the  bed  so  that  the  near  bed  turbulence  levels  are  reduced.  Clearly  such  a  difference  in  the 
development  of  turbulent  characteristics  of  the  boundary  layer  would  be  reflected  in  the  sediment 
mobilizing  and  transporting  capabilities  of  the  flow.  Such  a  difference  may  be  significant  even 
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for  fine  grained  sands.  Ii  is  also  important  to  note  that  the  results  presented  in  Conley  and  Inman 
[submitted]  were  for  a  smooth  and  rigid  permeable  bed.  It  is  not  clear  how  the  separate 
phenomena  identified  in  this  work  would  interact  and  reinforce  each  oilier  in  the  presence  of  an 
moveable  bed.  The  observation  of  an  asymmetry'  in  the  development  of  fluid-granular  boundary 
layers  discussed  in  Conley  and  Inman  [1992]  occurred  in  fine  grained  sediments.  It  was  shown 
that  this  asymmetry  was  not  due  to  asymmetries  in  fluid  velocity  or  acceleration  which  suggests 
that  the  combined  effects  of  ventilation  arc  indeed  significant,  even  in  fine  grained  sands. 

The  analysis  in  this  paper  of  the  relative  importance  of  the  net  stresses  due  to  wave 
asymmetry,  bottom  wind  and  ventilation  assumes  that  their  affects  arc  mutually  independent  of 
each  other.  This  assumption  is  questionable  as  all  affects  are  mutually  present  in  nature  and  may 
well  interact  in  ways  which  can  not  be  predicted.  As  an  example  of  a  likely  interaction,  consider 
that  Longuet-Higgins  [1958]  points  out  that  the  bottom  wind  results  from  alternate  boundary 
layer  thinning  under  the  crest  of  the  wave  and  thickening  under  the  trough  of  the  wave  The 
phase  between  the  vertical  velocity  associated  with  this  thickening  and  thinning  and  the  horizon¬ 
tal  velocity  in  the  boundary  layer  gives  rise  to  the  bottom  wind.  Conley  and  Inman  (submitted! 
have  shown  that  ventilation  gives  rise  to  a  similar  thinning  and  thickening  under  the  crest  and 
trough  respectively  and  it  is  likely  that  the  two  affects  may  be  mutually  reinforcing. 

5.  CONCLUSIONS 

A  field  experiment  designed  to  detect  the  presence  of  wave  driven  ventilation  in  perme¬ 
able  beds  was  performed.  Wave  driven  ventilation  was  tested  as  a  possible  mechanism  for  the 
maintenance  of  beach  slope  and  this  mechanism  was  compared  against  other  mechanisms. 

1.  Utilizing  Darcys  law,  ventilation  which  is  driven  by  progressive  surface  gravity  waves  is 


observed  in  the  permeable  beds  underneath  the  fluid  on  which  the  waves  propagate. 


1.  Progressive  surface  gravity  wave  driven  ventilation  in  the  permeable  bed  beneath  die 
fluid  on  which  the  waves  propogate  is  observed  through  the  utilization  of  Darcys  law. 

2.  The  magnitude  of  this  (low  is  well  predicted  by  linear  theory  given  that  the  local  effec¬ 
tive  bed  thickness  is  known.  The  effective  bed  thickness  is  assumed  to  relate  a  homogeneous  bed 
of  constant  permeability  to  the  the  true  local  stratigraphy. 

3.  The  net  stress  arising  from  ventilation  is  found  to  be  an  important  phenomena  for  shon 
period  waves  and  medium  to  coarse  grain  beds.  This  net  stress  can  dominate  the  stresses  arising 
from  velocity  asymmetries  in  shoaling  waves  or  steady  streaming  velocities  in  the  bottom  boun¬ 
dary  layer. 
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ences  under  grant  N00014-89-J-1060  with  the  University  of  California.  San  Diego.  The 
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Field  Experiment 


Schematic  of  the  experimental  layout.  Differential  pressure  sensor  is  buried  in  the  sand  with  pressure  port  separated 
in  the  vertical  only.  Absolute  pressure  sensor  and  em  current  meter  are  placed  in  the  water  column  directly  above 
the  differential  pressure  sensor. 
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Figure  5:  Solid  line  is  ventilation  spectra  as  calculated  from  differential  pressure  sen¬ 

sor.  Dotted  and  dashed  lines  represent  ventilation  spectra  predicted  from 
(11).  Dotted  line  assumes  bed  thickness  of  250  cm  and  dashed  assumes 
thickness  of  580  cm.  Units  are  in  (cm/s)2/cps. 
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Figure  6:  Schematic  illustrating  the  concept  of  effective  thickness.  Shallow  inhomogeneous  bed  of  sand  and  cobble  on  left 

can  be  modeled  as  homogeneous  bed  of  sand  on  right  with  effective  thickness  greater  than  the  true  thickness. 
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Figure  7:  Schematic  showing  the  effect  of  ventilation  on  oscillatory  boundary  layers. 

Flow  into  the  bed  enhances  bed  stress  while  flow  out  of  the  bed  decreases 
bed  stress.  In  an  otherwise  symmetrical  oscillatory  flow  the  result  is  a  bed 
stress  asymmetry  which  leads  to  a  net  stress.  The  direction  of  this  stress  is 
in  the  direction  of  the  orbital  velocities  experiencing  flow  into  the  bed. 
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Figure  8:  Comparisons  of  net  stress  arising  from  various  mechanisms  as  a  function  of 

mean  grain  diameter.  Solid  line  is  net  stress  due  to  ventilation  calculated 
from  (23).  Dashed  line  is  stress  due  to  bottom  wind  (21)  and  the  dotted  line 
is  stress  from  a  third  order  stokes  wave  (18).  All  values  are  normalized  by 
(pcfa2o2 /{2sinh2(kh)\).  8a  uses  wave  parameters  representative  of  field 
experiment  while  b.  and  c.  show  the  effect  of  longer  period  waves. 


74 


REFERENCES 

Batchelor,  G.  K.An  Introduction  to  Fluid  Mechanics,  615  pp.,  University  Press,  Cambridge, 
1985. 

Bijker,  E.  W.,  J.  P,  Th.  Kalkwijk,  and  T.  Pieters,  "Mass  transport  in  gravity  waves  on  a  sloping 
bottom",  Proc.  Conf.  Coastal  Eng.,  I4th(l),  447-465,  1974. 

Conley,  D.  C.,  and  D.  L.  Inman,  "Ventilated  oscillatory  boundary  layers",  J.  Fluid  Mech.,  submit¬ 
ted. 

Conley,  D.  C.,  and  D.  L.  Inman,  "Field  observations  of  the  fluid-granular  boundary  layer  under 
near-breaking  waves,"  J.  Geophys.  Res.,  97(C6),  9631-9643,  1992. 

Grant,  U.  S„  "Waves  as  a  sand  -  transporting  agent",  Amer.  J.  Science,  241(2),  1 17-123,  1943. 

Gu,  Z.,  and  H.  Wang,  "Gravity  waves  over  porous  bottoms".  Coastal  Engin.,  15, 497-524,  1991. 

Inman,  D.  L.,  and  R.  A.  Bagnold,  "Littoral  processes".  The  Earth  Beneath  the  Sea,  The  Sea,  vol  3, 
edited  by  M.  N.  Hill,  pp  529-553,  Interscience,  John  Wilet  &  Sons,  New  York,  1963. 

Koh,  R.  C.  Y.  and  B.  Le  Mehaute,  "Wave  shoaling,"  J.  Geophys.  Res,.  71(8),  2005-2012,  1966. 

Krumbein,  W.  C„  and  G.  D.  Monk,  "Permeablility  as  a  function  of  the  size  parameters  of  sedi¬ 
mentary  particles",  Am.  Inst.  Min.  and  Met.  Engin.  Tech  Pub  1492,  153-163,  1943. 

Longuet-Higgins,  M.  S.,  "The  mechanics  of  the  boundary-layer  near  the  bottom  in  a  progressive 
wave",  Proc.  Conf.  Coastal  Eng.,  6th,  184-193,  1958. 

Longuet-Higgins,  M.  S.,  "Mass  transport  in  water  waves”,  Phil.  Trans.  Royal  Soc.  London,  A, 
245,535-581,  1953. 

Putnam,  J.  A.,  "Loss  of  wave  energy  due  to  percolation  in  a  permeable  sea  bottom",  EOS,  Trans. 
AGU,  30,  349-356,  1949. 

Reid,  R.  O.  and  K.  Kajiura,  "On  the  damping  of  gavity  waves  over  a  permeable  sea  bed,"  EOS, 
Trans.  AGU,  38,  662-666,  1957. 


